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Abstract 


In  this  report  the  problem  of  robust  detection  of 
non-coherent  signals  in  noise  is  investigated  under  the 
assumption  that  the  noise  distribution  is  unknown,  but  , 

ji 

a member  of  a known  class  of  distribution  functions. 

This  problem  is  divided  into  three  main  categories; 
partially-coherent  signals,  unknown  frequency  and  phase, 
and  completely  random  phase. 

In  the  partially-coherent  problem,  two  general  methods 
are  introduced  for  the  design  of  robust  detectors  which 
guarantee  a non-trivial  lower  bound  on  the  probability  of 
detection  and  a non-trivial  upper  bound  on  the  probability 
of  false  alarm.  Two  recievers  are  designed  for  the 
special  case  in  which  the  distributions  of  the'  noise 
inphase  and  quadrature  phase  components  are  members  of  a 
class  of  p-point  distributions.  Simulation  results  for 
finite  sample  size  are  also  given  for  different  distributions. 

The  M-detector  method  is  extended  to  give  a general 
solution  for  the  second  problem. 

In  the  third  problem,  three  different  solutions  are 
presented.  All  consider  that  the  disttibution  of  the 
envelope  of  the  noise  is  a member  of  some  class  of 
distribution  functions  which  is  defined  by  quantiles. 

Simulation  results  for  these  three  detectors  are  given  at 
finite  seimple  sizes  for  different  distributions. 


:( 
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Throughout  the  report,  the  log  normal  and  contaminated 
normal  distributions  are  utilized  in  conducting  the 
simulations . 
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I . INTRODUCTION 

» 

I The  different  approaches  for  the  design  of  a receiver 

for  detection  of  signals  in  additive  noise  fall  between  two 
extremes.  On  one  hand  is  the  classical  parametric  approach 
which  requires  an  exact  and  complete  statistical  descrip- 
tion of  both  the  signal  and  the  channel  and  utilizes  the 
likelihood  ratio  as  a test  statistic.  On  the  other  extreme 
are  nonparcunetric  or  distribution  free  approaches  which 
require  minimal  information  about  the  utilized  channel  and 
lead  to  constant  false  alarm  detectors,  such  as  sign,  rank 
and  median  detectors  [1] . 

Often  in  practice,  the  information  available  is  only 
partial,  either  because  the  channel  description  is  incom- 
plete or  because  some  aspects  of  the  channel  vary  either 
temporally  or  spacially.  This  partial  information  is  usual- 
ly more  than  required  for  design  of  a nonparametric  detec- 
tor, but  is  insufficient  for  the  design  of  a parametric 
detector.  The  treatment  of  such  problems  depends  on  the 
nature  of  the  missing  information,  which  can  be  divided  into 
two  main  categories. 

1)  Parameter  Ambiguity:  In  this  case  there  is  a 

mathematical  expression  for  the  channel  model  with 
some  cimbiguity  about  the  exact  value  of  some  parame- 
, ters;  i.e.,  the  actual  noise  is  a member  of  a known  ij 

finite-dimensional  class  of  distribution  functions.  For 
example,  the  channel  is  normal  with  unknown  mean  and/or  | 

I variance.  ii 
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2)  Distribution  Ambiguity:  This  is  the  intrinsically 
more  difficult  case  where  there  is  no  mathematical 
expression  for  the  noise  distribution,  but  some 
information  about  the  distribution  is  available  such 
as  the  first  m moments.  This  means  that  the  noise 
distribution  is  a member  of  an  infinite-dimensional 
class  of  noise  distributions. 

The  problems  in  the  first  case  are  easier  to  handle 
than  those  in  the  second  and  in  some  particular  situations 
it  is  possible  to  achieve  an  optimal  performance.  At  other 
times  only  a minimax  or  adaptive  detector  is  possible. 
Actually,  the  problems  in  this  category  have  been  studied 
extensively  as  the  composite  hypothesis  problem  and  many 
good  adaptive  and  minimax  techniques  have  been  developed 
[2-4]. 

In  treating  the  problems  in  the  second  category,  the 
receiver  designers  divided  into  three  groups.  The  first  and 
most  conservative  group  accepted  the  nonparametric  tech- 
niques to  at  least  guarantee  the  probability  of  error  of 
the  first  kind  or  the  probability  of  false  alarm.  The 
second  group  utilized  the  available  partial  information  to 
obtain  the  best  fit  from  a selected  class  of  noise  distri- 
butions, then  designed  a parametric  detector  based  cn  this 
approximate  distribution.  Utilization  of  nonparcimetric 
techniques  in  such  cases  is  useful,  but  their  use  results 
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in  ignoring  most  of  the  available  information.  On  the  other 
hand,  generation  of  an  approximate  noise  distribution, 
although  desirable  since  it  leads  to  a parametric  receiver, 
may  be  very  dangerous  since  the  actual  observations  may  not 
come  from  the  assumed  class  of  distributions.  This  implies 
that  the  detector  performance  may  deteriorate  markedly  from 
the  expected  performance.  To  demonstrate  this  point,  let 
the  approximate  density  function  be  f ^ (x)  under  H ^ , then 

N 

the  likelihood  ratio  will  be  A(x)  = Jl  -x — , ^ . Now, 

i=l  0'’^i' 

fl(x  ) 

a single  factor  -e — j y equal  (or  almost  equal)  to  zero  or 

1 ' k^ 

infinity  will  upset  the  test  statistic  and  change  the  final 
decision  [5] . 

Out  of  these  two  groups  and  as  a result  of  observations 
similar  to  those  above,  came  the  third  group  which  rejected 
both  approaches  and  went  further  to  reject  the  utilization 
of  strictly  pareimetric  techniques  since  these  systems  are 
sensitive  to  bad  data  and  there  is  no  known  procedure  to 
protect  any  device  completely  against  the  presence  of  such 
outliers  or  bad  data  [6]  and  [7] . This  group  believe  the 
best  solution  to  the  problem  is  in  the  design  of.  detectors 
which  make  use  of  the  maximum  possible  amount  of  the 
available  information  and  perform  well  over  the  class  of  all 
distributions  which  possesses  the  properties  described  by 
this  information.  These  detectors,  which  are  called  robust 
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detectors,  are  not  expected  to  perform  as  well  as  the  optimal 
detectors  for  all  distributions  in  this  class;  however, 
they  should  perform  better  than  the  nonparametric  method 
for  a large  subset  of  this  class.  The  search  for  such 
robust  procedures  has  been  very  active  during  the  last 
few  years  [5,8-10],  All  of  the  above  work  assumes  that 
the  signal  is  coherent  and  all  except  [10]  assumes  some 
class  of  contaminated  distributions.  In  [10] , a general 
procedure  for  robust  detection  of  known  signals  in  additive 
noise  was  presented.  These  receivers  are  called  M-detectors. 

In  this  report,  the  problem  of  "Robust  Detection  of 
Non-Coherent  Signals  in  Additive  Noise"  is  considered.  In 
general,  the  problem  is  divided  into  three  basic  cases. 

1)  Partially-coherent  signals:  In  this  case,  it  is 
assumed  that  the  signal  has  an  unknown  phase  which 
is  constant  over  each  sequence  of  observations,  but 
changes  randomly  from  one  sequence  to  another. 

2)  Unknown  frequency  and  phase:  It  is  assumed  that  both 
the  signal  frequency  and  phase  are  unknown.  This  is 
usually  the  case  in  detection  in  the  presence  of 
Doppler  effects. 

3)  Random  Phase:  In  this  case  it  is  assumed  that  the 
unknown  phase  changes  randomly  from  one  sample  to 


another. 
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Chapters  II  and  III  investigate  the  first  problem. 

In  Chapter  II  we  extend  the  concept  of  M-detectors 
presented  in  [10]  to  the  noncoherent  case.  We  show  that 
under  some  sufficient  conditions  on  the  family  there 
exists  a transformation  which  maps  the  input  data  into 
some  intermediate  space.  Utilizing  these  points  aj  new 
input  data  we  construct  a detector  and  show  that  under 
some  assumptions  on  the  signal  amplitude  it  will  be 
robust  in  a max-min  sense  over  the  family  of  distributions 
under  consideration.  In  Chapter  III  we  present  another 
detector,  called  the  stochastic-approximation  detector, 
which  has  similar  properties  as  the  M-detector.  We  also 
give  simulation  results  to  show  the  finite  sample  size 
performance  of  the  above  detectors  and  to  compare  it  to 
the  conventional  square  law  detector.  Chapter  IV  extends 
the  results  of  the  previous  two  chapters  to  the  second 
case.  In  all  three  chapters  it  is  assumed  that  some 
information  about  the  inphase  and  quadrature  phase 
components  is  known. 

In  Chapter^  aset  of  detectors  is  presented  to  treat 
the  third  problem.  All  the  detectors  considered  in  this 
chapter  depend  on  the  envelope  of  the  received  observations 
as  input  data. 


In  Chapter  VI  the  performance  of  the  M-detector  with 
dependent  data  is  investigated. 
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Throughout  this  report,  the  stability  of  the  threshold 
with  respect  to  the  actual  noise  distribution  and  the 
power  of  the  test  were  used  as  criteria  for  judging 
receiver  performance.  Also,  for  all  simulations,  the 
noise  distribution  was  assumed  either  to  be  lognormal  or 
contaminated  normal.  These  two  distributions  are  the 
most  common  in  the  above  detection  problems  [11,121.  The 
normal  distribution  was  utilized  principally  for 
comparisons. 


I 
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II . 1 Introduction 

This  chapter  investigates  the  detection  of  signals 
with  unknown  phase  in  additive  noise  under  the  assumption 
that  available  information  about  the  noise  distribution 
is  not  complete.  Specifically,  we  shall  consider  these 
situations  where  there  is  no  available  mathematical 
expression  for  the  noise  distribution  function,  either 
because  it  changes  in  shape  with  uncontrolled  factors 
or  just  because  of  the  absence  of  enough  data.  On  the 
other  hand,  it  will  be  assumed  that  some  information 
is  available  about  the  distribution  of  the  inphase  and 
quadrature  phase  noise  components.  The  signal  component 
will  be  considered  as  samples  from  a demodulated  version 
of  the  received  signal  which  is  given  by 

Sj.(t)  = B Sin(w^t+0)  0<t<^T  (I.l) 

where  6 is  the  unknown  pahse  which  is  a random  variable 
with  unifoirm  distribution  over  the  interval  [-  n , IT  ] , 

T is  the  total  observation  period  andB  is  the  signal 
amplitude  which  is  equal  to  zero  under  and  is 

greater  than  zero  under  . We  shall  call  this  problem 

the  partially  coherent  case. 
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This  problem  appears  in  many  communication  systems 
when  there  is  uncertainity  about  the  jiiase  of  the  wave 
form  generated  by  the  transmitter's  ossilators.  In 
radar  applications,  it  appears  in  the  case  of  detection 
of  fixed  targets  [3  pp  335],  A more  important  aspect  of 
this  study  is  that  it  furnishes  a basis  for  the  more 
important  problems  in  radar  systems  such  as  the  design 
of  doppler  processors.  In  such  situations,  both  the 
frequency  and  phase  will  be  assumed  unknown. 

Usually  in  practice,  the  engineer  uses  some  version 
of  the  square  law  detector,  not  because  of  its 
optimality  but  because  of  its  simplicity.  In  this 
detection  procedure  the  detector  sums  the  inphase 
component  samples  and  the  quadrature  samples,  and  then 
utilizes  the  square  root  of  the  sume  of  the  squares  of 
the  above  two  sums  as  a test  statistic.  This  method 
has  two  main  defects.  First,  since  the  sum  is  a linear 
operation,  if  we  have  a large  noise  sample,  either  because 
the  actual  distribution  has  a long  tail  or  because  of 
some  measurement  error,  then  this  observation  is  more 
likely  to  upset  the  decision.  Second,  if  the  distribution 
is  not  known  exactly  then  setting  the  threshold  will 
require  an  adaptive  technique.  Usually,  these  adaptive 
techniques  cause  some  loss  in  the  power  of  the  test  [13]  . 
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In  the  following,  we  give  a procedure  for  the  design 
of  robust  detectors  that  guarantee  an  upper  bound  on 
both  probability  of  false  alarm  and  probability  of 
missing  under  the  above  mentioned  assumptions.  We 
design  one  of  these  detectors  for  the  class  of  p-point 
distributions,  where  only  one  quantile  is  known.  We 
compare  this  detector  to  the  square  law  detector  for  the 
lognormal,  conteiminated  normal  and  Guassian  distributions. 
The  lognormal  was  chosen  because  of  its  long  tail  and 
because  it  appears  in  many  practical  situations  such  as 
detection  in  sea  clutter  and  atmospheric  noise.  The 
contaminated  normal  is  shorter  in  tail  than  lognormal  and 
represents  sea  clutter  interference  in  some  sea  states 
[14]  . It  also  can  be  considered  as  a moderate 
representation  of  measurement  errors  [7]  . The 

Gaussian  was  chosen  because  the  square  law  detector  is 
optimal  against  it. 


II . 2 Problem  Statement 

Consider  the  set  of  observations  {z^}  given  by 
its  inphase  and  quadrature  phase  components  denoted  by 
{x^}  and  respectively  where 


Xi  = Ai  + w^ 


i = 1, . . . 


N 


Yi  = *2  - 


(II. 2) 
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where  {w^}  and  {v^}  are  two  sequences  of  independent 
identically  distributed  noise  which  are  also  mutually 
independent,  f(w^)  and  f(v^)  are  the  density  functions 


of  w^  and  respectively,  which  are  unknown  members 


of  a class  of  symmetric  density  functions  F , and  A2 
are  the  inphase  and  quadrature  phase  signal  components 


such  that  *^a2+a2  = a.  We  define  the  two  hypotheses 


Hq  and  as 


Hq:  a = 0 


(II. 3) 


A > 0 


We  are  interested  in  finding  a detection  procedure  to 


distinguish  between  and  such  that  both  the 


probabilities  of  error  of  the  first  and  second  types  have 
the  best  achievable  nontrivial  upper  bounds. 

In  the  following  analysis  we  shall  consider  f (v) 
and  f (w)  to  have  the  same  mathematical  expression  and 
we  shall  deal  only  with  the  case  where  A^^  and  A2  are 
constants  over  each  group  of  observations.  The  generali- 
zation to  the  case  where  A^^  and  A2  are  changing  in  a 
known  manner  as  functions  of  i is  straight  forward 
[10]  . 

Before  proceeding  to  the  solution,  we  shall  need  a 
few  additional  definitions.  Define  a class  of  functions 
C such  that  L e C if 


r 
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1)  L is  convex,  symmetric  about  the  origin  and  strictly 
increasing  for  positive  arguments. 

2)  l(t)  = dL(t)/dt  is  continuous  for  all  t . 

3)  for  all  f e F,  E^[l^(x)]  < “ 


4)  for  all  f e F,  9E^  [1  (x-9)]/9e 


exists  and  is 
non  zero  in  some  neighborhood  of  the  origin. 


Also,  define  as  the  value  of  6 which  minimizes 

NH 


N 


I L(x.-0) 
i=l  ^ 


and  as  the  value  of  6 which  minimizes 

NV 

N 

I L(y.-e) 

i=l  ^ 

or,  equivalently. 


(II. 4) 


(II. 5) 


(II. 6) 


and 


N 


I KVi-e^y)  . 0 
1 = 1 


and  define  0„  as 
N 


(II. 7) 


(II. 8) 


Define  also  B^(0|f)  and  B^(A|f)  as  the  probability  of 
false  alarm  and  the  probability  of  detection,  respectively, 
using  detection  strategy  d when  the  actual  noise  distribu- 
tion is  f ( • ) . 
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The  main  results  of  this  section  will  be  obtained  by 
verifying  the  following  two  inequalities; 

1 f e F (II. 9) 

and 

f e F (II. 10) 

where  B,  (vlf)  = lim  B, (Alf)  with  A = v//N  , and 
d d 

then  by  showing  that  BLQ(v|fQ)  is  related  to  the  power 
of  the  optimal  coherent  detector,  when  the  true  noise 
density  function  is  f^  , by  the  same  relation  as  that 
between  the  power  of  a coherent  and  noncoherent  detector 
when  the  noise  is  Gaussian.  The  subscript  Lq  refers 
to  a detector  based  on  a threshold  test  using  0^^ 
drived  from  e C as  a test  statistic;  i.e. , a test 

of  the  form 

“l 

/N  0-,  < Y (II. 11) 

N > 

Ho 

we  shall  call  this  detector  the  M-detector  on  components. 

In  fact  it  might  seem  more  appropriate  to  show  that 
BLp(v|fQ)  is  greater  than  or  equal  to  the  power  of  any 
randomized  test  of  hypotheses  under  the  above  assumptions 
and  assuming  that  the  true  noise  is  fQ(*)»  but  with  the 
absence  of  an  optimal  procedure  to  design  such  a detector,  it 
is  the  authors  opinion  that  the  proposed  criterion  is 
the  best  possible  one. 
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II . 3 Problem  Solution 

Assume  for  the  moment  that  the  class  of  density 
functions  F contains  a density  f (•)  of  minimum 
Fisher  information  for  location  , that  is,  a density 
for  which 

1(F)  = / [f  '(x)/f  (X)  f {x)dx 

— oo 

is  a minimum  over  all  f e F . Also,  assume  that  the 

* 

class  C of  functions  contains  a function  L such  that 

* * 

L (x)  = -log  f (x)  • 

* * 

For  this  case  we  will  show  that  (L  ,f  ) = (LQ,fQ);  i.e. , 
it  satisfies  equations  (9)  and  (10)  described  above. 

First,  we  require  a preliminary  lemma  whose  proof  can 
be  found  in  [15]  . It  is  also  a special  case  of 

Lemma  1 [10]  so  we  shall  only  outline  it  here. 

Lemma  ( 1 1 . 1 ) ; 

Whenevery  L £ C and  f e F , /n  (Oj^^-A^)  [/N  (ej^^-A2)  ] 
is  asymptotically  distributed  as  a zero-mean,  normal 
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Proof : 


1)  As  N A,  and  6^^  A,  a.s.  and  in 


NV 


probability. 


2)  Assume  that  A=0  then 

NH  - /jJ  /fj  1 


(11.13) 


N 


3)  P{—  y l(x. ^) } =>  Normal  (y,V^) 

/N  i=i  ^ ”•- 


N->-oo 


6=0 


as  Je  Ef  f 1 1 I 0= 

and 


= / l^(x)  f(x)dx 


From  the  above  lemma  and  under  the  ass\imption  that 
{w.}  and  {v.}  are  independent  of  each  other  we  obtain 


with  = t/N  (Gj^)  , 


, I . Tn  -2v2(f,L) 

fCT^lH  ) - 2 ~~  ® 

“ N->»  v^(f,L) 


and 


f(Tj^|Hi) 


llL 


2 2 
T +V 

" 2V^(f ,L) 


(11.14) 


) (11-15) 


V (f,L) 


v^(f,L) 


where  mbdefined  Bessel  function  of  the 

first  )cind. 

Then  the  asymptotic  power  of  a test  based  on 


satisfies 


B, (vjf)  = / 


pL_  e I,  (-4^— ) dT^ 

Y v^(f ,L)  ° v^(f,L)  ” 


= Q(oi,y) 


(11.16) 


where  Q{.,.)  is  Marcam's  Q function,  a = - ^ and 

v^(f,L) 

Y is  the  threshold  of  the  test.  Notice  that  for  any 
finite  falue  of  y » if  A is  a constant  then  v -*■  “> 
as  N -*•  “>  and  Q(a,Y)  1;  i.e.,  the  test  is  always 
consistent. 

In  the  above  discussion  it  was  assumed  that  and 

w^  are  independent  for  all  i.  Nevertheless,  the  results 
are  still  correct  if  v^  and  w^  are  given  by 


V.  = u.  cos  <p . 

11  1 

i = 1, . . . , N (11.17) 

w^  = u^  sin  <(>^ 


where  {u^}  and  are  sequences  of  independence 

identically  distributed  random  variables  with  (J)^ 
uniformly  distributed  between  [0,2n].  To  show  this, 
we  have 


P [e 


NH 


®NV 


= P{  i Z 1 
i=l 


N k . N k 

(X  - ^)<0,  i Z Ky  -^)<0> 
^ /N  /N  i=l  /N 


N 

Z x,<o. 


N 

Z Y,<0} 


1 


(11.18) 
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but  for  any  two  constants  a and  b 


[ — J X.  + — 'l  Y.]  is  asymptotically  distributed 

/N  i=l  ^ i=l  ^ 

1 ^ 

as  a normal  random  variable,  then  [ — I X.]  and 

/n  i=l  ^ 

1 ^ 

[ — y . ] are  asymptotically  jointly  normal  random 

/N  i=l  ^ 

variables  [16]  and  it  is  enough  to  show  that  they  are 
uncorrelated  to  prove  their  independent.  Since 

E{(  i I X.) ( i I Y.)} 
i=l  ^ /n  i=l 


1 N 

= ^ I E{X.Y.)  + (N-1)  E(X.)E(YJ  (11.19) 

^'-—1  JLi  X X 


then 


cov  [ { i I X . ) , ( - I Y . ) ] 

i=l  /n  i=l 


1 N 

^ I E(X.Y.)  - E(X.)E(Y.) 


- »<s> 


(11.20) 


and  since  the  variance  of  each  is  of  order  1 then  the 
correlation  coefficient  is  of  order  (^)  and  goes  to  zero 
as  N ^ i.e.,  they  are  asymptotically  independent. 


Then 


17 


N k 

=>*■'  J H-i-  ;^>i<» 


p{i  ,I  i(yi-4>lo 

/N  i=l  ^ 


= p(/N  6^^<k^)  P(/I 


(11.21) 


and  the  results  of  the  previous  discussion  continue  to 


hold.  This  completes  the  proof  of  the  following  lemma 


Lemma  I I . 2 : 


Whenever  L e C and  f e F,  v^  independent  from 


for  all  i=l,...,  N or  given  by  eqn.  (11.17), 


will  be  asymptotically  distributed  under  as 

m.  2 


f(T  ) = 2 ■ e 

^ v^(f,L) 


■"N 

2v2(f ,L) 


(11.22) 


and  under  H, 


(Tjj2+v2) 

2v2(f,L) 


f(T  ) = -2  - e 

^ v^(f,L) 


T V 

I ( N ) 

o'  2,. 

V (f,L) 


(11.23) 


To  study  the  power  relations  for  different  f e F 


we  need  the  following  lemma. 


Lemma  (II . 3)  : 


If  there  exists  an  f^  e F such  that  LQ=-log  f^  e C 


then 


sup  Bj^(v|fQ)  = BLQ(v|fQ) 
Lc  G 


(II. 24) 


J 
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and  if 


v^(f ,LQ)<v^(fQ,LQ)  for  all  f e F 


(11.25) 


then  the  following  relation  holds  asymptotically 


BLo(0|f)  < BLQ(0|fQ)  for  all  f e F 


(11.26) 


Proof : 


From  equation  (11.22) 


Blo«>|£)  = '2v2(f,L„)l 


< exp  i 


-v2 


2^(fO/Lo)  ^ 


by  virtue  of  equation  (11.25). 


Also,  for  any  level  of  false  alarm  a and  for  any  L e C, 
if  fg  is  the  true  density 


= -2v^(f-,L)  In  a 7^2 

0 


(11.27) 


T„  V 


“ T 2v^(f0,L) 


Y V (fQ»L) 


O'  2,.  "N 


V-2  Ina 


oo  1 9 V 

/ ' V(fo,L) 


R e 


In(R 


V 


0'*'  v(fQ,L) 


)d  R 


= Q[ 


v(fQ,L) 


, /-2  Ina  ] 


(11.28) 


but  as  v^(f«,L)  decreases,  v,  = ^ v increases  and  so 


1 V(fg,L) 


J 
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L 


2 2 

does  Bj^(v|fQ).  On  the  other  hand  v (fQ,LQ)£v  (f^/L) 


since  this  is  the  variance  of  the  maximum  likelihood 
estimates  of  and  A2  and  is  equivalent  to  the 

Cramer-Rao  lower  bound.  This  completes  the  proof  of 
the  lemma. 


To  get  the  relation  between  BLQ(v|f)  and  BLQ(vlfQ) 
we  assume  that  BLQ(0|fQ)  = a,  then 


BLn(v|f)= 


I 


v(fQ,LQ)/-2  Ina  v (f/LQ) 


2 ^ 2 

Z Z + V 1 

— — = ] 


2v‘(f ,Lq) 


zv 


0 2,_  _ . 

V 


) dz 


using  R = 


v(f ,Lq) 


BLn^vlf)  = 


/ Re 

v(ff.,L^)  

Y '■= — rs— 7 — 5 — Ina 

v(f ,L-) 


R^v*2 


Iq(Rv-)  dR 


= Q(v' ,y') 


(11.29) 


where  v'  = 


v(f  ,Lq) 


and 


Y-  = I ^ 

Wf^T^ 


Notice  that  as  v{.,.)  decreases  both  v'  and  y'  increases 
with  the  Scime  ratio,  which  makes  this  relation  untractable 
analytically,  especially  with  the  absence  of  an  analytic 


J 
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expression  for  the  Q function  itself.  On  the  other  hand, 
we  can  write  eq  (11.29)  in  the  form 


Q(v' ,y' ) = Q(r  v,r  y) 


(11.30) 


v(fo'^0^  - V 

^ = v(f7L~)-  i V = , and  if  V > Y 


then  [ 3 ] 

2 2 

“(v*  "^"Y*  ) °° 

Q(r  v,rY)  = 1 - e I (— ) "^  I (V  y' ) 

n=l  v‘ 


2-2 

-r  V /I  .,,2. 

— ^ (1+k  ) “ - 

= 1 - e ^ I (k)  I (rv^k) (11.31) 

n=l  " 


Y 

where  k = — < 1.  Taking  the  derivative  with  respect  to  r. 


^ 2 

^ Q(r  v,r  Y)  = G[Ij^(r^  v k)-k^  lQ(r^  v 'k)  ] (11.32) 


where  G is  a positive  constant.  Thus 

? 2 

3 - I,  (r  V k)  2 

Q(r  v,r  Y)  > 0 iff  ^ • This  implies 

d JT  — t V ““ 

lQ(r  v k) 

that  if  V and  y satisfy  the  above  inequality  for  all 
values  of  r^l  then  the  power  of  the  test  will  be 
monotonic  in  r 


1 BLo^vIfg) 


(11.33) 


This  completes  the  proof  of  the  following  lemma. 
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Lemma  (II . 4)  ; 

If  there  exists  fQEF  such  that  Lo=-log  igCC  and 

v^(fo'LQ)  > v^(f,LQ)  f e F (11.34) 

I^(r^  v^k)  2 _ 

and  if  k for  all  r>l  and  v>y  then 

lQ(r^  v^k) 

f E F.  (11.35) 

The  above  inequality  was  solved  for  values  of  y^I.9226 

which  is  equivalent  to  Pp£.15  and  it  was  found  that  k 

should  be  less  than  or  equal  to  (0.91166).  Values  of 

V are  such  that  the  probability  of  detection  is  greater 

than  or  equal  to  about  60%.  Notice  that  this  is  the 

region  of  interest  in  most  communication  applications. 

Notice  also  that  in  the  case  of  robust  coherent  detectors 

[10  ] this  range  was  for  all  Pj^.  5 . 

The  above  lemma  gives  the  lower  bound  on  the 

probability  of  detection  over  the  class  F.  To  investigate 

how  good  this  lower  bound  is,  we  compare  the  proposed 

M-detector  to  the  asymptotically  most  powerful  coherent 

detector  when  the  true  noise  density  is  fQ(.).  If  we 

assume  that  the  probability  of  false  alarm  is  a , for 

the  asymptotically  most  powerful  coherent  detector 

1 

Bopt(Vilfo)  = 1-  {$'^(l-a)-Vj^  (11.36) 


( 
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where  4*  ( . ) is  the  standard  normal  distribution  function 

and  I(.)  is  the  Fisher  information  number  for  location. 

2 1 

For  the  proposed  M-detector,  since  v (f ^ ,Lq)  then 


- z^+V2  Kfo^ 


/-2  Ina 


= Q(v2  I (fg)  ,/-2  Ina) 


Iq[z  V2  I (fg)! 


(11.37) 


Comparing  the  above  two  equations  with  the  case  when 
f^(.)  is  Gaussian  [17]  we  find  that  if  both  probability 
of  false  alarm  and  detection  are  constant  then  — is 
the  same  in  both  cases,  i.e.  in  replacing  a coherent 
scheme  with  an  M-detector,  the  loss  in  signal  to  noise 
ratio  will  be  the  same  in  the  M-detector  case  as  in  the 
Gaussian  case. 

All  the  above  can  be  concluded  in  the  following 
theorem. 

Theorem  II. 1; 

If  there  exists  q(*)  e F 3 LQ=-log  f^  e C and 

v^(fg,LQ)  > v^(f,LQ)  Vf  e F 

then  the  M-detector  will  satisfy  the  following  relations 
asymptotically 

BLo(0|f)  lBLQ(0|fQ)  Vf  e F (11.38) 

Bl  (v|f)  < BLQ(v|fQ)  f L e C (11.39) 


% 

J 


and  if  v = 
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threshold  and 


Ij^(r^v  Y) 
lo^r^v  Y) 


> (-!-) 


Y,2 


for  all  r^l 
then 

BLo^vlf)  i BLQ(v|fQ)  ¥f  e F (11.40) 

The  above  results  show  that  the  design  of  the 
receiver  will  depend  on  finding  the  least  favorable 
density  function  f^  as  defined  in  Lemma  II. 3.  To  find 
this  density,  we  notice  that  it  must  have  a finite 
Fisher  information  number  for  location  according  to 
assumption  3 on  C.  Thus,  it  will  also  be  least  favorable 
in  the  family  G of  density  functions  defined  as 

G = {f:  f e F and  1(f)  < (11.41) 

However,  the  least  favorable  density  from  this  class  is 
the  one  with  minimum  Fisher  information  (cf.  Theorem  2 [15]). 
Thus,  to  find  f^,  we  follow  the  steps  described  below. 

a)  Find  f such  that  I(f  ) £ 1(f),  all  f e G.  This  can 

usually  be  done  using  Lagrange  multiplier  techniques. 

* 

b)  Check  that  -log  f e C. 

* 2 * 

c)  Check  that  f maximizes  v (L  ,f)  over  F.  In  general, 

* 2 * 

if  f does  not  maximize  v (L  ,f)  for  all  f e F,  then 


♦ 


f^  will  not  exist  since  in  this  case  we  have  only  two 
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alternatives.  Either  there  is  no  f which  maximizes 
the  variance  over  F,  or  there  exists  an  f which 
maximizes  the  variance  but  is  not  in  G,  so  L=-log  f 
is  not  an  element  of  C. 

In  summary,  the  results  given  above  show  that  the 
most  robust  detector  may  be  designed  by  finding  f^ 
maximizing 

v^(f,-log  f)  = (I(f))"^,  all  f e F,  (11.42) 

and  then  using  LQ=-log(fQ)  to  define  ^0  ^ 

basing  the  threshod  detector  on  0^. 

II . 4 The  Family  of  P-Point  Distributions 

This  section  deals  with  the  above  detection  problem 
when  the  noise  has  a density  f which  is  a member  of  the 
fcimily  F of  p-point  distributions  defined  as 

F = {f:  f(x)d(x)  = p, 

~ a 

f symmetric  and  continuous  at  +ci}  . (11.43) 

Reasons  for  choosing  this  specific  family  of 
distributions  are  as  follows. 

1)  This  family  covers  a very  wide  class  of  distributions; 
for  fixed  a and  p,  F contains  a scaled  version  of 
almost  every  symmetric  distribution. 
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2)  The  required  information  for  this  family. 


f(x)dx,  is  one  of  the  most  easily  measured 


parameters  of  a distribution. 


To  apply  our  procedure  to  this  family,  we  must  find  f. 


and  Lq  which  are  given  by  the  following  lemma. 


Lemma  (II. 5) 


Over  the  family  F of  density  functions  given  by 


(11.43),  the  density  function 


bj^  cos^  (Cj^x)  , 


X < a 


-cJx| 


(11.44) 


X > a. 


which,  for  appropriate  choices  of  bj^,b2»Cj^,  and  C2f  is 


continuous,  and  has  a continuous  first  derivative,  has 


the  following  properties:  a)  = -log  f^  e C,  b)  inf^^^^ 


v2(fQ,L)  = v2(fQ,LQ)  = l/I(fQ)=sup^^pv2(f,LQ) . The 


proof  of  this  lemma  can  be  found  in  either  [1£3,[13  or  in 


Lemma  II. 5 shows  that  the  detection  strategy  will 


be  to  find  0„„  and  such  that 

NH  NV 


I l(x^-6jjj^)  0 

1=1 


(11.45) 


I KyrV  = “ 
1=1 


(11.46) 
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Thus , we  have 


®Lo(vlf)  i fiLpCvlfg)  = Q(^  ,/-2  Ina) 


(11.54) 


II . 5 Numerical  Results 

To  exaunine  the  finite  sample  size  performance  of 
these  M-detectors,  we  apply  the  proposed  procedure  to 
3 different  distributions 
1)  Lognormal  distribution 


f(u)  = 


~ 2 2,1/2 

(2Tra  u ) ' 


exp 


where  o =6  db. 

2)  Contaminated  normal  distribution 

2 2 

f(u)  = d-r"^)  exp  (— ~^— ) + o 

2o  g 


-u 


exp( — 
2g 


2r(l-r)  ,-u  (g  +1),  ,u  (g  -l)v 

— exp  (- 22 ' ^0 

go  4 g a 


2,  2 
t (g  - 
T . 
4 g a 


Notice  that  this  is  the  probability  density  of  the 
envelope  of  two  orthogonal  components,  each  having  a 
density  function  of  the  form 

(1-r) 


2 _ 2 

f (w)  = ^ — 5 T*7T  2 2 1/2  ^ 2 2 ^ 

(27ro'^)-^'^  2o^  (2Tig"^o"^)  2 g^o^ 


r and  g was  taken  to  be  0.25  and  2.25  respectively. 
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3)  Rayleigh  distribution 

2 

f(u)  = u e ^ 

Each  of  the  above  distributions  was  considered  to  be 
a member  of  a p-point  class  of  distributions.  We  set 
p=.5,  then  found  a.  Having  a and  p the  pareimeters 
of  the  nonlinearity  lj^(.)  in  equation  (11.50)  were 
calculated,  then  equations  (11.48)  and  (11.49)  were 
utilized  to  find  6„  . Using  a test  of  the  form 
given  by  (11.51)  was  conducted.  To  find  the  value  of 
the  threshold  y a run  was  made  first  under  using 

10^  groups  of  samples  then  a minimum  mean  square  error 
technique  was  applied  to  fit  the  output  data  to  a curve 
of  the  form 

- In  Pp  = ao  + ai  Y + a2  Y^. 

The  above  equation  was  then  used  to  set  the  threshold 
for  low  probabilities  of  false  alarm. 

Figures  1 and  2 show  the  probability  of  detection 
Pj^  verses  the  signal  to  noise  ratio  (S/N)  for  different 
values  of  probability  of  false  alarm  For  log-normal# 

A 

(S/N)  here  is  defined  as  20  log^^Q  j;j,  where  M is 

the  median  of  the  distribution  of  the  envelope. 

Figures  3 and  4 show  the  same  for  the  contaminated  normal 
and  the  Rayleigh  distributions.  Notice  that  the  signal 
to  noise  ratio  required  to  achieve  a certain  probability 


S/N  (dB) 


probability  of  detection  vs  S/N  for  the 
log-normal  distribution  (o=6  dB)  and  M-detector 
on  components:  N=30. 


Figure  1 


S/N  (dB) 

probability  of  detection  vs  S/N  for  the 
log-normal  distribution  (o=6  dB)  and  M-detector 
on  components:  N=10. 


Figure  2 


S/N  (dB) 

probability  of  detection  vs  S/N  for  the 
contaminated  normal  distribution  ( r=0. 25,g=2. 25) 
and  the  M-detector  on  components;  N=10 


Figure  3 


S/N  (dB) 

probability  of  detection  vs  S/N  for  the 
Rayleigh  distribution  ^ the  M-detector  on 
components;  N=10. 

Figure  4 
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F 


of  detection  at  any  fixed  value  of  probability  of  false 
alarm  in  the  case  of  lognormal  is  higher  than  that  for 
contaminated  normal,  which  in  its  turn  is  higher  than 
the  S/N  in  the  Rayleigh  case.  This  is  due  to  the  fact 
that  the  lognormal  has  the  longest  tail  and  the  Rayleigh 
has  the  shortest  tail. 

Figures  (5-8)  show  verses  S/N  for  a square  law 
detector  for  the  above  distributions.  Figures  (9-12) 
compare  the  performance  of  the  conventional  square 
law-detector  to  the  performance  of  the  proposed  M-detector. 
Examining  figures  (9)  and  (10)  we  notice  that  the 
M-detector  always  has  a better  performance  than  the 
square  law  detector.  Notice  also  that  the  improvement 
in  S/N  becomes  better  as  P decreases  and  achieves  about 

r 

4 db  at  10  For  the  contaminated  normal,  figure  (11) 

shows  that  this  improvement  is  between  (.5)  and  1 db 
with  the  M-detector,  better  than  the  square  law  detector. 

It  should  be  taken  into  account  that  the  p-point  family 
is  not  the  best  fcimily  to  represent  the  contaminated 
normal  and  much  more  improvement  can  be  achieved  by  the 
utilization  of  some  other  distribution  classes  e.g.  21. 

For  the  Rayleigh  distribution,  figure  (12)  shows  that 
the  M-detector  is  about  1 db  worse  than  the  square  law 
detector.  This  loss  should  be  expected  since  the  square  law 
detector  is  optimal  in  this  case. 


S/N  (dB) 


probability  of  detection  vs  S/N  for  the  contaminated 
normal  distribution  (r=0 . 25 ,g=2 . 25)  and 
the  square  law  detector;  N=10 

Figure  7 


S/N  (dB) 


Comparison  of  square  law  detector  ( ) and  M-detector 

on  components  ( — ) in  log-normal  (*r:=6dB)  interference;  N=30 

Figure  9 


2 


6 
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S/N  (dB) 

Comparison  of  square  law  detector  (- 
on  components  ( — ) in  log-normal  (P=6dB) 

Figure  10 


-)  and  M-detector 
interference;  N=l0 


-2  0 2 4 6 8 

S/N  (dB) 

Comparison  of  square  law  detector  ( ) and  M-detector 

on  components  ( — ) in  Rayleigh  interference;  N=10 


Figure  12 
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From  the  above  results  we  can  conclude  that  utilization 
of  minimal  information  such  as  one  quantile  might  lead  to 
improvement  in  the  performance  of  about  4 db  or  more. 
Comparing  this  to  a loss  of  1 db  in  the  Rayleigh  case,  and 
taking  into  account  that  the  longer  tail  distributions  are 
much  more  likely  than  the  shorter  tail  ones,  we  find  the 
M-detector  is  highly  preferable  than  the  square  law 
detector.  One  should  also  expect  to  get  more  improvement 
in  the  performance  of  the  M-detector  by  utilization  of 
more  information.  For  example,  if  we  have  some  idea 
about  some  of  the  partial  moments,  the  family  of  p-point 
distributions  would  be  replaced  by  the  generalized  moment 
constrained  class  defined  by 

F = {f:  /S„  fdx  < P , l<n<N 

or 

F = {f:  fs  fdx  = P , l<n<M  and 
■'  n n — — 

fs  fdx  < P , M+l<n<N} 

Notice  that,  if  S (x)  = 1 for-a<x<a  and  zero  otherwise 

n — 

we  get  a quantiles  and  if  ^ -b£x£b  zero 

otherwise  we  get  partial  moments.  More  details  about 
this  class  are  in  [22] . 

An  important  remark  about  these  similation  is  that 
if  we  normalized  the  data  with  the  appropriate  value  of  a. 


A 


\ 

I 


then  we  can  set  one  threshold  to  guarantee  upper  bound 
on  the  probability  of  false  alarm  even  at  small  values 
of  N as  N=10.  Also,  a lower  bound  on  the  probability 
of  detection  exists.  Ihese  bounds  usually  are  the  thres- 
hold and  signal  required  to  achieve  certain  Pp  and 
when  the  distribution  is  lognormal.  In  fact,  it  was 

noticed  that  the  changes  in  P„  and  P_  when  the 

r D 

distribution  of  the  noise  changes  from  lognormal  to 
contaminated  normal  are  very  small. 

II . 6 Summary 

In  this  chapter  we  presented  a procedure  of  the 
detection  of  signals  with  unknown  phase  in  noise,  under 
the  assumption  that  the  noise  distribution  is  unknown 
member  of  a class  of  symmetric  density  functions  F.  It 
was  shown  that  this  method  asymptotically  guarantees 
an  upper  bound  on  both  probability  of  false  alarm  and 
missing  whenever  the  noise  density  function  is  element 
in  the  class  F under  consideration.  Simulation  results 
were  given  to  show  the  finite  sample  size  performance 
of  this  method  and  to  compare  it  the  conventional  square 
law  detector. 
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III . Detection  of  Partially  Coherent 
Signals  in  Noise  II 


111.1  Introduction 

In  the  previous  chapter  the  problem  of  detection  of 

partially  coherent  signals  in  noise  was  solved  by 

introducing  the  M-detector  on  components.  In  this 

chapter,  we  continue  investigating  the  same  problem  as 

defined  in  Section  II. 2.  We  introduce  another  approach 

which  is  called  the  stochastic  approximation  detector 

(SA) . It  may  be  noticed,  that  the  M-detector  depends 

on  finding  6„„  and  0„„  as  M-estimates  of  A,  and 
NH  NV  1 

A_  respectively,  and  then  conducting  a threshold  test 

2~ 

using  0„  = 0„„+0„,,  as  a test  statistic.  The 
N NH  NV 

SA-detector  will  be  exactly  the  same  as  the  M-detector 

except  that  the  M-estimates  (0„„  and  0»_J  will  be  replaced 

NH  NV 

by  their  stochastic  approximation  counterparts  S„„  and 

Nrl 

Sjj^  respectively. 

In  the  following  section  we  shall  give  a brief 
description  of  the  SA-estimates  and  their  properties. 

111. 2 Stochastic  Approximation  estimation 

Consider  the  sequence  of  observations  {x^}  given  by 

Xi  = s + Wi  i = 1,...,  N (III-l) 


1 

J 


I 
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where  {w^}  is  a sequence  of  independent  identically 
distributed  random  variables  with  a density  function 
f(w),  which  is  an  unknown  member  of  a class  of 
symmetric  density  functions  F,  and  S is  an  additive 
signal.  It  is  of  interest  to  find  an  estimate  S„  of 
the  form 


= S 


N-1 


N ^ (N-1) ^ 


(III-2) 


where  B is  a constant  and  g{*)  is  some  function 

such  that  S„  is  a consistant  estimate  of  S and 
N 

2 

asymptotically  normal  with  variance  v (f,g).  It 
would  be  of  interest  also  to  find  a pair  (f^jg^)  such 
that 

Vg(fo»go)  “ Min  Max  Vg(f,g)  (III-3) 

g f eF 

This  problem  has  been  investigated  by  several 
researchers  and  we  shall  summarize  some  of  their  results 
here. 

Assume  that  y (d)  = E[g(x-d)]=0  has  a unique 
solution  d=s.  Consider  also,  the  following  set  of 
assumptions ; 

Assumption  (A^)  : y is  a Borel-measu'rable  function;  and 
(d-s)y(d)  > 0 


for  all  d 5^  s. 
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Assumption  (A2) : For  some  positive  constants 
and  for  all  d 

k|d-s|<|y(d)  |_<k^|d-s| 


k and  k,  , ;| 

i-l 


Assumption  (A^) : For  all  d 
M (d)  = n (d-s)  + 6 (d, s) 

where  6(d,s)  = 0(|d-sl)  as  (d-s)  -*■  0 and  n -*■  0 


2 

Assumption  (A.):  a)  sup  E [g  (x-d) -y  (d)  ] <«> 

^ d 

b)  lim  E [g (x-d) -y (d) ] ^=0^ 
d^s 

Assumption  (A^) : 

2 

lim  lim  sup  / z (d)dR=0 

R-y°°  ld-s|<e  |z(d)|>R 

where  z (d)  = g(x-d)  - y (d) 

Theorem  (Ill-l)  [23]; 

Suppose  that  assumptions  (A^ ) through  (A^)  are 

V2  ^ 

satisfied  and  2Bk>l.  Then  N '^(S^-S)  is  asymptotically 

normally  distributed  with  mean  0 and  variance 

2 2 -1 
B o (2B  n -1) 

The  proof  will  be  omitted.  Sacks  also  showed  that 
the  above  theorem  remains  true  (with  K replaced  by  n) 
if  assumption  (A2)  is  replaced  by  (A^'): 
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Assumption  (A2'):  For  ail  d and  some  positive 
1 y :d)  l<kj^ld-s  j 

and  for  every  that  0<tj^<t2<“ 

inf  1 M (d)  I >0 

til|d-s|<t2 

The  above  theorem  gives  conditions  under  which  the 
SA-estimate  is  consistent  and  asymptotically  normal. 

To  achieve  the  minimax  property,  define  the  class  of 

* 

density  functions  F as 

F*={feF:  I(f)<“}  . (III-4) 

* 

Assume  that  there  exists  a density  f_(.)eF  such  that 

^ f 

I(f,.)<Iff)  for  all  f in  F and  define  ( . ) =-:^( . ) / then 

0 - U Iq 

we  have  the  following  theorem  from  [22] 

Theorem  (III-2)  ; 

Suppose  that  F is  a convex  set  of  probability 

* 

density  functions  with  F non-empty,  gQ(-)  satisfies  the 
conditions  of  theorem  (III-l)  and 


m'f{a)  = 3^  E [-g^  (x-a)  ] =-/  <3qM  f'(x+a)dx  (III-5) 

Then 

Vs(f,go)lv2(fo,go)<v2(fQ,g  ) (111-6) 

for  all  f(.)  e F 


48 


Notice  that  the  above  conditions  are  very  similar  to 
the  required  conditions  for  the  M-estimates.  To  see  the 
difference  in  performance  between  the  two  estimators  we 
use  the  following  lemma  from  [22]  . 


Lemma  (III-l)  ; 

If  the  M-estimate  and  SA-estimate  are  both  consistent 
and  asymptotically  normal  at  a density  function  f with 
positive  finite  variance  then 


v^(f,LQ)  < Vg(f,gQ) 


(III-7) 


with  equality  iff  m^ (a) =Ef ^ [g^ (x-a) ] . 


The  above  means  that  the  M-estimate  will  always  have 
a lower  or  equal  variance.  Another  problem  here  is  Sq, 
the  initial  value  of  to  be  used  in  equation  (III-l)  . 
Although  will  not  depend  much  on  the  specific  value 

of  Sq  asymptotically,  this  choice  will  effect 
dramatically  at  finite  sample  sizes.  Notice  also  that  if 

f=fQ  and  g=go=-f^  ' = TTf^  if  B = , 

then  ~ ^ which  is  the  Cramer- Rao  lower 

bound  and  the  variance  of  the  estimate  of  S. 

In  the  following  section,  we  apply  the  above  results 
to  the  detection  problem  in  Section  II. 
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III. 3 Stochastic  Approximation  Detectors 

Consider  the  problem  of  detection  of  partially 
coherent  signals  described  in  Section  (II. 2)  where  {x^} 
and  inphase  and  quadrature  components  of 

the  observations,  respectively,  given  by 


X.  = A,  + w. 
Ill 


i=l, . . . , N 


(III-8) 


Yi  = *2  Vi 


with  the  same  assumptions  and  definitions  in  Section  (II-2) 


Define  S 


and 


Sj^  as  the  stochastic  approximation 


estimates  of  A^  and  A2  respectively  given  by 


H 


Si  = S(i_i)  + 


B 


H 


f 9("i  - 


(III-9) 


and 


Si  = ^Ll)  T^^^i  - ^I-l)) 


(III-IO) 


Define  S . as 
1 


s..  = 

1 


(iii-ii) 


and  the  SA-test  of  hypotheses  as 

Hi 

< y 


where  y is  the  threshold  of  the  test. 


. (III-12) 
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where  B^(0|f)  [Bg(v|f)]  refers  to  the  probability  of  false 

alarm  [detection]  using  the  SA-detector  with  nonlinearity 
g when  the  actual  noise  distribution  is  f . 

The  above  theorem  shows  that  whenever  the  M-detector 
on  components  guarantees  upper  bounds  on  the  two  kinds  of 
error  probability,  the  SA-detector  will  guarantee  the 
same  bounds.  On  the  other  hand,  for  any  density  function 
for  which  the  performance  is  better  than  these  bounds 
the  M-detector  will  have  better  performance.  This  means 
that  in  using  the  SA-detector  in  place  of  the  M-detector 
we  should  expect  a slight  deterioration  in  performance 
in  return  for  the  computational  convenience.  In  the 
next  section  we  shall  give  some  simulation  results  to 
evaluate  the  finite  sample  size  performance  of  the 
SA-detector  and  to  find  out  how  much  loss  in  efficiency 
we  suffer  by  replacing  the  M-detector  by  the  SA-detector. 

III-4  Numerical  Results 

To  evaluate  the  finite  sample  size  performance  of 
the  SA-detector,  simulations  were  done  for  three  different 
distributions,  lognormal,  contciminated  normal  and  Rayleigh. 
These  distributions  are  the  same  utilized  in  the  study 
of  the  M-detector.  To  conduct  these  simulations,  each 
distribution  was  considered  to  be  a member  of  a class  of 
p-point  distributions  with  p=0.5  as  before.  Then  fQ(-) 
is  given  by 


fo(x)  = 


cos  (Cj^x)  |x|  < a 


-C2X 


(III-17) 


go(x)  = 


b2  e 


b tan{Cj^x) 


X > a 


X < a 


(III-18) 


b tan  (Cj^a)  sgn  (x)  I X I > a 


with  b^,b2/b,Cj^  and  C2  chosen  such  that  is  in  F 
and  gQ(»)  is  continuous  at  +a.  The  constant  B in 

equation  (III-l)  was  chosen  to  be  equal  to  ^ . 

2 ® 

As  a result,  Vg(f,gQ)  should  be  asymptotically  equal 

to  B for  all  f in  F.  Three  different  choices  for  Sq, 

the  initial  value  of  the  estimate,  were  utilized. 

1 \ qH  — yr  qV  = u 


1)  = Xj^  and  ^ ^1  • 

2)  Sg  = Sj  = 0 . 

N 

I X. 

.•Ill  1 


3)  S 


H _ i=l 


N 

A.  ^i 


and  S, 


The  third  choice  led  to  the  best  performance,  so  we 
report  here  on  the  results  of  these  simulations  only. 

Figures  13,14,15  and  16  show  the  probability  of 
detection  for  the  SA-detector  against  the  signal  to  noise 
ratio  (S/N)  for  lognormal  N=30,  lognormal  N=10,  contaminated 
normal  N=10  and  Rayleigh  N=10,  respectively.  Figures  17, 

18  and  19  compare  the  performance  of  the  SA-detector  to 
the  M-detector  for  lognormal  and  contaminated  normal. 


S/N  (dB) 


probability  of  detection  vs  S/N  for  contaminated 
normal  distribution  (r=0. 25,g=2. 25)  and  SA-detector;  N=10 

Figure  15 


S/N  (dB) 


Comparison  of  SA-detector  ( — 
in  log-normal  (a 


-)  and  M-detector  on  components  ( — ) 
=6  dB)  interference;  N=30 

Figure  17 


A 


Comparison  of  SA-detector  ( ) and  M-detector  on  components  ( — ) 

in  log-normal  (o=6  dB)  interference;  N=10 


Figure  18 
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and  show  a loss  of  up  to  about  2db  in  signal  to  noise  ratio. 
Figures  20,21,22  compare  the  SA-detector  to  the  square  law 
detector.  Although  the  SA-detector  is  better  for  the 
lognormal  case  it  looses  its  superiority  in  the  case  of 
contaminated  normal,  and  runs  about  Idb  worse  than  the 
square  law  detector. 

III-5  General  Discussion: 

The  above  simulations  show  that  the  deterioration  in 
performance  as  a result  of  replacing  the  M-detector  by 
the  SA-detector  is  not  slight  but  serious.  It  also  shown 
that  the  performance  of  the  SA-detector  is  affected  not 
only  by  the  choice  of  Sq,  but  also  by  the  first  few 
observations.  To  clarify  this,  consider  N=10.  Then  if 
the  first  two  observations  are  large  with  the  Scime  sign, 
then  the  chance  of  correcting  this  effect  by  the  other 
eight  observations  is  zero.  On  the  other  hand,  the 
M-detector  can  be  well-approximated  by  a one  step 
Newton-Raphson  version  [20] . In  addition,  the 
M-detector  may  also  be  approximated  by  an  on-line  system 
if  we  are  not  interested  in  the  specific  value  of  0^^ 
but  only  in  its  relation  to  the  threshold  y . By  the 
monotonocity  of  the  nonlinearity  1(.)  in  its  argument, 
if  than  either 


S/N  (dB) 


Comparison  of  SA-detector  ( — ) and  square  law 
detector  { ) in  log-normal  interference  (o=6  dB) ;N=30 


Figure  20 


2 


4 


6 


8 


S/N  (dB) 

comparison  of  SA-detector  ( — ) and  square  law 
detector  ( ) in  log-normal  interference  (o=6  dB) ; 

Figure  21 


2 


6 


8 


S/N  (dB) 

Comparison  of  SA  detector  ( — ) and  square 
in  contaminated  normal  interference  (r= 

Figure  22 


l^w  detector  ( 
.25,g=2.25)  ; N 
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N 

y l(x.-d)  > 0 
i=l 

or  (III-19) 

N 

y K-x.-d)  > 0 
i=l 

and  if  6„>y  then  there  exists  a number  d such  that 
N 

■“  *“2  2 2 

3nd  |0,_,)>d,  where  d'^=Y  -d  . To  construct  the 
' HN ‘ — ' VN ‘ 

above  approximate  detector,  we  pick  a set  of  numbers 
{Yji^}  between  y//2  and  y and  conduct  the  test  as 
in  equation  (III-19)  . Note  as  the  number  of  {Yj^}"*"® 
we  get  the  exact  M-detector.  An  example  of  the  above 
implementation  is  shown  in  figure  (23) . 

III. 6 Summary 

In  this  chapter  we  analysed  another  method  for 
detection  of  signal  with  unknown  phase  in  additive  noise 
called  stochastic  approximation  detector.  It  was  shown 
that,  although  this  method  is  very  attractive  and  has 
good  properties  from  the  theoretical  point  of  view,  it 
fails  completely  in  application  at  finite  sample  sizes. 

In  addition,  an  approximate  on  line  implementation  for 
the  M-detector  on  component  was  presented. 
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IV.  Detection  of  Signal 
with  Unknown  Phase  and  Frequency 

IV. 1 Introduction 

In  the  previous  two  chapters,  the  problem  of  detection 
of  signals  with  unknown  phase  in  the  presence  of  additive 
noise  was  considered  under  the  assumption  that  the  noise 
distribution  function  is  not  completely  known.  It  was 
assumed  that  all  other  signal  parameters  — except  perhaps 
the  signal  amplitude  — are  known.  Often  in  practice,  in 
addition  to  the  phase,  the  signal  frequency  is  not 
completely  known.  The  analog  received  signal  in  such 
cases  can  be  modeled  as 

s^(t)  = B(t)  sin[(w^+w)t  + 4)]  (IV. 1) 

where  w^  is  a known  carrier  frequency.  4)  is  an 

unknown  phase  taking  any  value  in  [0,211]  and  w is  an 
unknown  frequency  in  the  interval  [w',w"].  Although 
this  problem  appears  in  many  communication  applications 
as  a result  of  the  uncertainty  about  the  transmitter 
oscillator  frequency,  it  is  more  frequently  encountered 
in  radar  applications,  especially  in  detection  of  moving 
targets  [24].  In  radar  applications,  this  problem  is 
usually  referred  to  a^' Doppler  detection"  and  w is  called 
the  "Doppler  Shift".  Since  our  concern  in  this  chapter 
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will  be  mainly  with  the  radar  problem,  it  is  in  order  here 
to  discuss  briefly  the  nature  of  this  Doppler  shift. 

Assume  that  there  exists  a target  which  is  moving 
with  radial  velocity  v , then  the  range  of  this  target 
as  a function  of  time  will  be  given  by 

R(t)  = Rq  - V t (IV-2) 

where  Rq  is  some  constant  equal  to  the  initial  range  at 
t=0.  If  the  transmitted  signal  is 

s^(t)  = B(t)  cos  (w^t)  , (IV-3) 

then  the  signal  returned  from  the  target  will  be 

Sj.(t)=  B{t-  T(t))  cos  [w^(t-  T(t))]  (IV-4) 

where  t (t)  is  the  round-trip  delay  time 
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0 

where  x = — ^ and  is  constant.  Equation  (IV-7)  shows 
that  the  velocity  of  the  moving  target  will  have  two 
main  effects  on  the  reflected  signal; 

a)  A shift  of  the  carrier  frequency. 

b)  A compression  or  stretching  of  the  time  scale  of  the 
eve lope. 

The  first  effect  is  usually  more  serious  than  the 
second  one  [24] . In  our  study  here  we  shall  not  take 
the  second  effect  into  consideration  since  the  generaliza- 
tion of  the  solution  will  be  straightforward. 

Before  proceeding  to  the  solution  it  would  be  helpful 
to  discuss  the  case  in  which  the  additive  noise  is  known 
to  be  Gaussian. 

IV. 2 The  Gaussian  Noise  Case 

Consider  the  envelope  of  the  transmitted  signal  to  be 
a sequence  of  pulses  each  with  duration  D;  the  separation 
between  pulses  is  T.  Then  if  the  pulse  amplitude  is  A and 
the  observations  are  corrupted  by  additive  Gaussian  noise, 
these  observations  will  be  in  the  form, under 

r^  = A cosl(w^+w)t  + 4)]  + n(t) 

Hq.*  r^  = n(t) 


i=0, . . . , (N-1)  (IV-8) 


i 


But  since  the  noise  is  Gaussian  and  independent,  then  the 


conditioned  likelihood  ratio  would  be  expressed  as 


(N-1) 

A (r|  w,((>)  = n X (v.  |w  , (()) 

^ _ A 1 


(IV-9) 


where 


X (r^|w,(J)) 


-1  2 
exp{  — ^ / [r.(t)-A  cos((w  +w)t+4))]  dt} 

2a  iT  ^ ^ 


iT+D 


exp{  / [r^(t)  1 dt} 

2a  iT 


iT+D 


= exp{  — 2 / r^(t)  A cos  [ (w  +w)  t + (|)]dt 

a^  iT  ^ 


iT+D 


, j.  J.T  2 2 

j / A COS  [ (w  +w)t+(|)]dt  (IV-10) 

2o  iT  ^ 


Substituting  from  (IV-9)  in  (IV-10) 


N-1  iT+D 


X(r|w,(|))  = F exp{  — =r  I f r.(t)A  cos  [ (w  +w)t+(})]dt 

2a  i=0  iT  ^ ^ 


(IV-11) 


where  F corresponds  to  the  received  energy  and  is  not  a 


function  of  w or  4).  Next,  we  break  the  remaining  cosine 


terms  up  according  to  the  trigonometric  identity 


cos[WQt+(j))  = cos  (w^t)  cos  (4))  - sin(wQt)  sin  (4)) 


J 


1 


r 
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iT+D 

/ r.(t)  cos [ (w  +w) t] dt  (IV-12) 

1 c 

iT+D 

/ r. (t)  sin[(w  +w)t]dt  (IV-13) 

iT  1 c 

Thus,  we  can  write  equation  (IV-11)  as 

A(x|w,({i)  = F exp{  — ^ [L  cos  ip-L  sin  (fi]  (IV-14) 

_ ^ c s 

0 

If  we  assumed,  as  is  usual  in  radar  detector,  that  (fi 
is  uniformly  distributed  over  the  interval  [0,211],  then 

A (r|w)  = [A  (r  |w,(j))  ] 

2tt 

= f Tir  F exp{  —j  [L  cos  4i-L  sin  4i]  } d([) 

0 ZTT  c s 

1/2 

= F Iq[A{L^+lJ)  /o^]  (IV-15) 

where  modified  Bessel  function  of  order  0. 

Notice  that  equation  {IV-15)  constitutes  a complete  solution 
to  the  detection  problem  when  w is  known  [3].  If  the 
frequency  is  unknown,  as  is  the  case  here,  we  face  one  of 
two  situations.  Either  we  are  interested  in  the  specific 
value  of  w or  not.  If  we  are  not  interested  in  the  value  of 
w,  we  can  average  equation  (IV-15)  with  respect  to  w to 
get  A(r)  and  then  compare  it  to  a threshold  [17].  If 


and  define 


N-1 

Lc=  I 

^ i=0 


and 


N-1 

^s  = ^ 

® i=0 
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we  are  interested  in  w,  then  we  choose  a set  of  possible 


values  for  w,  say  , 

a fine  mesh  on  the  interval 


such  that  this  sequence  forms 
tw',w"].  Next,  we  formulate 


a multiple  hypotheses  problem  with  a hypothesis  for  each 
of  the  M-values  of  w . That  is, 

j 


r^  = n(t) 


r^  = A cos  [ (w^+w^^)  t + (})]  + n(t) 


H.. : r.  = A cos  [ (w  +w„)  t + ({)]+  n(t) 

Ml  c M 

i = 0,. . . , (N-1)  (IV-16) 

The  likelihood  ratio  comparing  the  iOj  hypothesis 
to  the  null  hypothesis  is 

= X(xlw^)  i = 1,...,  M (IV-17) 

If  no  X^  is  greater  than  the  threshold,  then  the  null 
hypothesis  is  accepted;  otherwise,  each  w^  for  which  the 
threshold  is  exceeded  is  a potential  detection.  If  only 
one  target  is  present,  the  maximum  likelihood  estimate  of 
w will  be  that  w^  for  which  X^  is  the  largest. 

In  the  above  analysis,  it  was  assumed  that  the  noise 
distribution  is  Gaussian  with  known  parameters.  If  the 
variance  of  this  Gaussian  distribution  is  unknown,  then  an 


V.- 


r 


71 


adaptive  technique  will  be  necessary  to  set  the  threshold. 
Many  types  of  such  techniques  are  possible  [25] . One  of 
these  methods  is  the  one-parameter  threshold.  In  this 
procedure  we  utilize  a set  of  reference  cells,  say  K cells, 
then  we  integrate  the  pulses  in  each  of  these  cells  to 
get  z.,  j=l,...,  K.  The  threshold  is  then,  with  q constant, 
k 

T = q I ^ . (IV-18) 

j=l  K 

Although  this  method  is  simple  enough,  it  has  two  main 
disadvantages:  a)  it  depends  heavily  on  the  assumption 
of  normality  which  is  doubtful  in  most  radar  application, 
b)  if  there  is  another  target  in  one  of  the  reference 
cells,  this  will  increase  the  threshold  and  original  target 
might  be  missed. 


IV. 3 M-poppler  detector 

Examining  the  previous  detection  scheme,  we  notice 
that  it  consists  mainly  of  two  processes.  First,  assuming 


that  the  true  frequency  is  w^  , we  transform  the  original 


set  of  observations  r^(t)  into  a new  set  of  discrete 


in-phase  and  quadrature  phase  components  x^ ^ and  y^^ , 
respectively,  given  by 


X . 


iT+D 


iT 


r.(t)  A cos  [ (w  +w.)t]dt 
1 c 3 


(IV-19) 


and 
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iT+D 


= / r- . (t)  A sin  [ (w  +w . ) t]  dt 

-<  iT  ^ c ] 


iT 

Then,  under 

iT+D 


X.  • = / n(t)  A cos  ( (w  +w.)t]dt 

J-J  iT  c j 


(IV-20) 


w.  . 
ID 


(IV- 21) 


and 


iT+D 


IT 


n(t)  A sin  [(w  +w.)t]dt 
*-  D 


= V.  . . 

ID 


Under  H, 


iT+D 

Xi  • = / A cos  [ (w  +w)  t+i))]  A cos  [ (w  +w.)t]dt+w.. 

-*iT 


(IV-22) 


2 iT+D 


/ cos  [ (w  -w. ) t+())]dt  + w.  . (IV-23) 


iT 


c D 


ID 


and  since  w<<  — , then 
D 


X.  . ^ cos  [ (w-w  . ) iT+(|)]  + w.  . 

+D  ^ D ID 

Likewise , 


(IV- 24) 


y.  ■ sin[  (w-w. ) iT+4)]  + v.  . 

J-J  ^ J ID 


(IV-25) 


Noting  that  if  w=w.  the  first  term  in  both  (IV-24)  and 

^2  2 

(IV-25)  reduces  to  — ^ cos  <t>  and  — ^ sin  <|)  , respectively. 
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and  that  w. . and  v. . are  both  Gaussian  since  n(t)  is 
ID  ID 

Gaussian  by  assumption,  the  second  step  in  the  previous 
scheme  was  to  utilize  a square  law  detector  to  decide 
whether  a signal  is  present  at  w ^ . 

As  was  discussed  in  Chapter  II  of  this  report,  the 
performance  of  the  square  law  detector  depends  on  the 
assumption  of  normality  of  n(t),  and  if  n(t)  is  not 
normal,  as  is  the  case  of  interest  here,  the  performance 
of  the  square  law  detector  might  deteriorate  seriously. 

To  avoid  this,  we  propose  a more  robust  doppler  detector 
which  is  called  the  M-Doppler  detector.  This  detection 
scheme  consists  also  of  two  steps.  The  first  is  exactly 
the  same  as  in  the  previous  scheme.  The  second  utilizes 
the  sequences  ^^ij^  inputs  to  a set  of 

M-detectors  on  components  presented  in  Chapter  II.  Since 
the  mathematical  work  required  for  this  procedure  is  the 
same  as  in  Chapter  II,  we  are  going  to  give  only  a brief 
description  of  the  procedure. 

a)  Find  the  nonlinearity  1(.)  which  is  optimum  in  the 
minimax  sense  over  the  family  of  distributions  F 
which  contains  f(w^j)  and  f(v^j). 

b)  Find  e„t,(j)  and  (j)  defined  as 

NH  NV 


N-1 


Jo  ‘ “ 


(IV-26) 


and 
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N-1 

I 

i=0 

llVij  - 

(IV-27) 

Find 

d)  Compare  0j^(j)'s  to  a threshold  y 

It  should  be  noticed  here  that  an  advantage  of  this 
method  is  the  stability  of  the  threshold  with  respect  to 
the  distribution  of  n(t),  so  an  adaptive  technique  is 
not  needed  and  there  is  no  consequent  loss  in  the  power 
of  the  test.  Figure  (24)  shows  an  M-Doppler  detector. 


IV- 4 Numerical  results 

To  test  the  performance  of  the  proposed  detector 
a set  of  simulations  was  conducted  under  the  assumption 


that  w. . 

1] 


and  v^j  are  components  of  a lognormal 


distribution  with  a = 6 db  , and  then  repeated  for  the 

case  when  they  are  contaminated  normal  with  r = .25  and 

g = 2.25.  In  both  cases,  f(w^j)  and  f(Vj^j)  were  considered 

to  be  members  of  a class  of  p-point  distributions  with 

p=0.5  . Then  1(.)  was  found  to  be  of  the  form 
r tan (ct)  I t I £ a 

Kt)  = ) 

{^tan  (ca)  sgn  ( t)  |t|  > a 


for  appropriate  c such  that  l(t)  is  a continuous  function. 

The  set  w^  was  chosen  such  that  the  interval  [w' ,w"] 
is  divided  into  16  equal  parts.  Using  the  above  nonlinearity 


and  values  of  w ^ , the  simulations  were  conducted 


and  the 
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r^(t) 


i V'  ” 

C)^  ■ 

estimator 

ACOS  [ (W^+Wj^)t] 


[ 

vn 

M 

C)^ 

J 

estimatoi 

A sin  [ (W^+Wj^)  t] 


A sin[  (W(,+Wj^)t] 

S/N  (dB) 

M-Doppler  detector 
Figure  24 


Choose  the  largest 


76 


results  are  suii'jnarized  in  figures  (25)  and  (26)  which  show 
the  probability  of  detection  verses  the  signal  to  noise 
ratio  for  different  values  of  probability  of  false  alarm 
for  N=10.  An  important  result  which  is  not  shown  in  these 
graphs  is  that  the  simulation  confirmed  again  the  stability 
of  the  threshold.  It  was  found  that  if  we  normalized 
the  thresholds  in  both  cases  by  the  value  of  a,  we  get 
almost  the  same  results. 

Figures  (27)  and  (28)  compare  the  M-Doppler  detector 
to  the  one-parameter  adaptive  detector  discussed  in 
section  (IV-2) . It  may  be  noticed  that  the  superiority 
of  the  M-Doppler  detector  increases  as  P_  decreases  and/or 

r 

increases  and  the  difference  in  (S/N)  comes  to  about 
8 db  at  P_=10  ^ and  P =.99  for  the  lognormal  case,  and  more 
than  this  for  the  contaminated  normal  case. 


S/N  (dB) 

Comparison  of  M-Doppler  detector  ( — ) and  one-parameter 
adaptive  detector  ( ) for  log-normal  distribution  (a=6dB) ;N=10 

Figure  27 
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! 

i 


S/N  (dB) 


[ 


Comparison  of  M-Doppler  detector  ( — ) and  one-parameter 

adaptive  detector  ( ) for 

contaminated  normal  distribution  {r=0. 25,g=2. 25) ;N=10 


t. 


Figure  28 
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V.  Detection  of  Signals  with  random 
Phase  in  Sea  Clutter 

V-1  Introduction; 

In  the  previous  part  of  this  report  we  assumed  that 
the  signal  phase  is  an  unknown  nonrandom  parameter  over 
each  set  of  observations  but  can  change  randomly  from 
one  set  to  another.  In  this  chapter,  we  move  one  step 
ahead  to  assume  that  the  signal  phase  changes  randomly 
from  one  observation  to  another.  The  main  difference 
between  these  two  cases  can  be  seen  if  we  consider  the 
distribution  of  the  inphase  and  quadrature  observation 
components.  In  the  first  case,  the  inphase  and/or  the 
quadrature  components  will  have  the  same  distribution 
under  both  hypotheses  except  for  a location  shift 
which  is  constant  or  almost  constant  under  all  observa- 
tions, so  it  was  suitable  to  utilize  a test  which  was 
sensitive  to  location  changes  to  distinguish  between 
the  two  hypotheses.  In  the  second  case,  the  distribution 
of  these  components  will  also  be  the  same  except  for 
a location  shift,  but  this  location  shift  will  change 
randomly  from  one  observation  to  another  and  will  have 
a mean  which  is  equal  to  zero  under  both  hypotheses.  This 
means  that  utilization  of  a hypothesis  test  which  is  sensitive 
! to  the  means  of  the  components  should  not  be  expected  to 

t 

lead  to  good  performance,  especially  if  this  test  is 


L 


the  envelope  of  the  observations  is  monotonically 
increasing  with  the  signal  in  the  case  of  random  phase 
[26  Figures  (4-7)].  As  a result  of  the  above 

observation,  the  main  trend  in  treating  the  problem  of 
detection  of  signals  with  random  phase  was  generally  to 
design  detectors  which  use  the  envelope  of  the  observa- 
tions as  input  data  instead  of  the  components  or  the 
observations  [17,  27  and  28  ] . _ In  addition,  if  we 


examine  the  results  in  the  previous  chapters  of  this 
report  we  may  notice  that,  for  the  types  of  noise  under 
considerations,  censoring  the  observations  always  leads 
to  improvement  in  the  performance.  This  conclusion  was 
also  supported  by  some  previous  results  by  Trunk  [27] 
when  he  tested  the  trimmed-mean  detector. 


As  a result  of  the  above  discussion,  we  are  going 
to  consider  in  this  chapter  only  these  detectors  which 
utilize  the  envelope  of  the  observations  as  input  data. 
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4 

1 i 

1i 

I 


and  which  limit  the  effect  of  the  tails  of  the  observations  ^ 

distribution  in  some  way.  It  should  be  emphasized  here  • 

that  we  are  looking  for  detectors  which  have  good  performance  | 

over  a wide  class  of  noise  distributions  when  there  is 
only  partial  information  about  these  distributions.  Also, 
we  shall  consider  that  the  available  information  is 
always  quantiles. 

In  Section  (V-2)  we  present  the  "M-detector  on 
envelope".  In  Sections  (V-3)  and  (V-4)  we  discuss  the 
"Robust  Quantizer"  detectors  and  the  "Extreme  Value 
Theorem"  detectors  respectively.  Section  (V-5)  gives 
brief  description  of  the  "Trimmed  Means"  detector  which 
[ will  be  utilized  in  Section  (V-6)  for  comparisons. 

r 

( (V-2)  M-Detector  on  Envelope: 

[ “ ■ — - 

I Let  Zj^  be  independent  identically 

' distributed  samples  taken  from  the  probability  density 

I f(z)  and  let  F(z)  be  the  distribution  function  of  z. 

[ 

I Consider  also  that 

/ 2~  2 , . 
z.  = x.+y.  (V-1) 

1 1 ■'i 

where 

Xj^  = A cos  + w^  (V-2) 

j and 

y^  = A sin  4)^  + v^  (V-3) 

I 

i 
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where  and  are  noise  components,  ij)^  is  a random 
phase  with  uniform  density  function  over  the  interval 
[0,2-n],  and  A is  the  signal  amplitude  which  is  equal  to 
zero  under  and  greater  than  zero  under  . The 
exact  value  of  A may  be  unknown.  Assume  also  that 
there  exist  two  known  real  numbers  and  Z2  such 

that  0<Z^<Z2<“  and  F(Z^|Hq)=Pj^  and  F(Z2|Hq)=P2  where 
Pj^  and  P2  are  known  and  not  equal  to  zero  or  one.  It 
is  required  to  design  a detector  which  gives  good 
performance  against  all  distributions  having  the  above 
properties,  and  which  is  insensitive  to  the  variations 
of  the  actual  distribution.  As  was  discussed  in  the 
previous  section  this  detector  should  be  expected  to 
limit  the  effect  of  large  observations. 

Before  proceeding  to  describe  the  M~detector  on 
envelope,  we  need  the  following  definitions.  Define 
the  nonlinearity  1(.)  as 


Kt)  = 

where  a 
before. 


r-h 


(t)  <z 


V. 


and  b are  two  constants  and  Z 


Define  6„  as  the  value  of  6 
N 


<t<Z2  (V-4) 

<t 

■^>^2  defined 

such  that 
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N 

I g(z  -0  )=0  (v-5) 

i=l  ^ 

g(z-0)  = l{z-0)  - E[l(z)  (Hq]  (V-6) 

Then  the  M-detector  on  envelope  will  be  of  the  form 
>1 

0j^  < Y (V-7) 

where  y is  the  threshold. 

Examining  the  above  detection  strategy,  we  find 

that  the  main  idea  is  to  use  a fraction  of  the  small 

observations  (those  less  than  Z^)  to  get  rid  of  the 

large  observations  (those  larger  than  Z_)  without  sorting 

^2 

the  data.  If  b was  chosen  to  be  equal  to  , then 

1 

asymptotically  this  will  be  equivalent  to  the  average 
sum  of  all  observations  in  some  intermediate  interval 
of  length  (Z2~Z^).  In  this,  the  M-detector  shares  with 
the  median  detector  the  good  property  of  deleting 
outlyers;  however,  it  does  not  depend  on  one  observation. 

So  it  is  midway  between  the  mean  and  median  detectors. 

To  study  the  stability  of  the  threshold  we  must 
first  obtain  an  expression  for  the  distribution  of  0^ 
under  Hq. 

Distribution  of  0„ 

Lemma  (V.l) 

Under  Hq,  ^ ^ 
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Proof : 


For  any  positive  number  c the  law  of  large  numbers 


implies 


1 N 

lim  - I g(z^-e)  = E[l(z  -e)|HQ]  - E[l(z  )1Hq] 


Nh-co  ” i=l 


since  l(x)  is  monotonic  in  its  argument.  Also 
1 ^ 

lim  - I g(z.+e)  > 0 

N->-oo  i=l 


But  since  6^^  is  defined  such  that  ^ g(z^-9^)=0,  then 

i=l 

for  any  positive  number  e , there  esits  a number  Nq  such 
that  for  all  N_Nq 


P{  I ©N  I— 1 ^ 


Lemma  (V. 2) : 

Under  Hq,  /N  0^  is  asymptotically  distributed  as 
normal  with  zero  mean  and  variance 


^2( El3_UU ^ 

[E'  (9(2-6) 


(V-8) 


Proof : 


^ ) < 01 


(V-9) 
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Using  the  central  limit  theorem  it  can  be  shown  that 
the  sum  in  equation  (V-9)  tends  to  normal.  Also, 

lim  E[g2(z  - ^ - E[g2(z)] 

^-►00  ^ /N 


lim  /N  E[g(z-  — )]  = lim  /n  E[1{z-  — ) 


- l(z)] 


--Kff 


= -k  E’  lg(Z-0) |g,„l 


Then, 


k E'  [g(z-0)  I ] 

P{/n  e.<k}  <!>  { 

E ' [g  (z)  ] 


(V-10) 


where  $(.)  is  the  standard  normal (distribution  function. 

To  evaluate  the  sensitivity  of  the  threshold  with 
respect  to  the  distribution,  we  calculate  the  explicit 
value  of  the  variance  in  equation  (V-10) . 


E'{g(z-e) = E'{l(z-e) 


(b+l)a 


f f ( z)  dz 
Z, 


= rp  _p  ] ■ 


(V-11) 


i.e.  it  is  independent  of  the  specific  shape  of  f(z) 
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E[g  (z)]  = 


E[l^(z)  ] - E^[l(z)  ] 


= [l-P2]  + [ba+^^i^  Z ^ [P^-P^]  [ l+P^-P^] 


+ [ 


-2 


{b+l)a,  ? 


Z 9 Z 1 _ 

2 1 

(b+l)a 


2^2  9 ^2 

It/  zf(z)dz-[J  zf(z)dz]'^} 


ni.p,-p,i  / zf(.,a. 

2 12  1 7 , 


(V-12) 


It  can  be  seen  that  the  value  of  the  threshold  will 
change  with  changes  in  the  true  distribution  only  through 
the  last  two  terms  in  equation  (V-12).  On  the  other  hand, 
since  we  are  dealing  with  detection  in  sea  clutter,  where 
all  possible  noise  distributions  range  between  Rayleigh 
and  lognormal  [29]  a proper  choice  of  Z and 

will  minimize  the  effect  of  these  two  terms.  For  example, 
if  we  choose  and  Z^  such  that  P^  and  P2 equal  .1  and 

.5  respectively,  the  effect  of  these  terms  on  the 
variance  will  be  negligible.  This  choice  also  would 
help  in  calculating  E[1(z)|Hq]  which  is  required  to 
perform  the  test. 


Receiver  implementation 

To  perform  the  above  test  as  it  was  described  in 
equation  {V-7) , we  must  first  calculate  0^.  This 
process  may  not  be  easy  if  the  niimber  of  observations 
involved  is  large  sine  1(.)  is  a nonlinear  function. 
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g{.)  is  monotonic 

such  that  £ > 0„ 
N 

N 

I g(z.-e)  < 0 
i=l 

and  if  e < e„,  then 
N 

N 

I g(z.-e)  > 0 
i=l 


But,  on  the  other  hand,  noting  that 
in  its  argument,  for  any  number  e 


Accordingly  if  we  determined  the  threshold  to  b^-  y , 
then  the  M-detector  on  envelope  can  be  implemented  in 
the  form 


N 


I 

i=l 


g 


< 


H, 


0 


(V-13) 


or,  equivalently 


I l(z.-Y)  < E[l(z) |h  ] (V-14) 

i=l  Hq 

The  implementation  of  (V-14)  is  shown  in  figure  (28*). 


Numerical  results; 

To  examine  the  finite  sample  size  performance  of  the 

M-detector  on  envelope,  simulation  experiments  were 

conducted  at  sample  sizes  of  N=10  and  30  for  the  case 

when  the  noise  has  a lognormal  distribution  given  by 

In^  z 

£(^,  , ^ 

/2ttz 


«D>A079  47S  JOHNS  HOPKINS  UNIV  BALTIMORE  MO 

robust  noncoherent  detection. (U) 
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B 


I 


where  o = 6 db,  and  at  N=10  for  the  contaminated  normal 


distribution  given  by 


-z 


9 5 

2 2 r^z^  2k' 

f{z)  = (1-r^)  e e 


^ 2_a-r)z 

^ 4k  ^ 4k 


(V-16) 


where  Iq(*)  is  the  modified  Bessel  function,  r=0.25 
and  k=2.25. 


In  each  experiment,  the  values  of  and  were 
found  such  that  P,  and  equal  0.1  and  0.5  respectively. 

1 p ^ 

* 2 

b was  taken  to  be  ^ = 5.  To  calculate  E[l(z) |Hq],  a 

Rayleigh  approximation  of  the  form 

- (z-m) ^ 

f(z)  = e (V-17) 


was  utilized.  First,  we  calculate  m and  a such  that  this 
distribution  satisfies  the  requirements  at  Zj^  and  Z 2-  Then, 
we  use  this  approximation  to  calculate  E[1(z)|Hq].  It 
was  found  that  the  maximum  difference  between  the 
approximate  and  true  values  is  about  5%.  It  was  noticed 
also  that  the  value  of  m is  very  small  in  both  cases. 

Figures  (29-31)  give  the  results  of  these  simulations. 
The  signal  to  noise  ratio  was  taken  as  20  log^^^  (A/M) 
where  A is  the  signal  amplitude  and  M is  the  median  of 
the  distribution  under  consideration.  Notice  that  the 


S/N  (dB) 

Probability  of  detection  vs  S/N  (dB)  for  contaminated 
normal  distribution  (r=0 . 25 ,g=2 . 25)  and  M-detector  on 
. envelope;  N=10 

Figure  31 


required  (S/N)  to  achieve  a certain  at  the  same  P„  in 
the  case  of  lognormal  is  slightly  higher  than  that 
required  in  the  contaminated  normal  case.  The  following 
remarks  can  be  made : 

^'l 

1)  — is  almost  the  same  for  the  two  distributions. 

2)  At  N=10,  if  the  thresholds  were  normalized  by  the 

appropriate  value  of  then  the  thresholds  which 

— 2 —6 

would  give  a P^,  of  10  and  10  when  the  noise  is 

-2  -7 

lognormal,  would  give  a P_  of  10  and  10  when 

r 

the  noise  is  contaminated  normal.  This  confirms 
the  stability  of  the  threshold. 

3)  At  N=10,  using  the  same  normalization  we  found  that 
at  the  same  signal  to  noise  ratio  the  probability 
of  detection  is  also  almost  constant,  independent 
of  the  true  distribution. 

Remarks  2 and  3 imply  that  the  detector  is  robust 
or  insensitive  to  changes  in  distribution,  even  at  small 
sample  sizes  such  as  N=10.  More  about  this  detector 
will  be  presented  later  in  this  chapter. 

Summary 

In  this  section,  the  M-detector  on  envelope  was 
presented  as  a solution  to  the  problem  of  detection  of 
signals  with  random  phase  in  sea  clutter.  The  main 
idea  behind  this  detector  was  the  use  of  some  of  the 
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very  small  observations  to  delete  the  large  ones  without 
need  for  sorting  the  data.  The  robust  properties  of 
this  detector  were  shown  both  theoretically  and 
experimentally.  In  addition,  a simple  method  for 
implementing  this  detector  on  line  was  presented. 

(V-3)  Robust  Quantizer  Detectors; 

The  main  idea  behind  the  set  of  detectors  to  be 

discussed  in  this  section  is  in  the  utilization  of  a 

quantizer  to  transfer  the  input  data  to  a set  of 

preassigned  values,  then  to  utilize  the  quantizer  outputs 

as  an  input  to  the  decision  rule.  Assume  that 

N m 

{z.}  is  the  sequence  of  input  data  and  {a.} 

^ i=l  ^ j=0 

is  a sequence  of  real  number  such  that  aQ=0  and 

a^=oo  . Then  the  output  of  the  quantizer  q^  will  be 

given  by 


q.  = q(z.)  = 1.  a.  ,<z.<a.  (V-18) 

j j“l“  1 j 

m 

where  {1.}  is  the  set  of  preassigned  values  for  the 

^ j=l 

quantizer  output.  The  test  of  hypotheses  will  be  of 
the  form 


H, 


1 N “1 

Tn  = ^ I < y 

i=l  ^ Hq 

where  y is  the  threshold.  The  above  choice  of 


(V-19) 


N 


in 


fact  does  not  in  anyway  limit  the  generality  of  the 


I 


procedure  since  the  q^^'s  are  independent  as  long  as  the 


z.'s  are,  and  we  are  still  free  to  choose  the  l^'s  to 
1 D 

satisfy  any  optimality  criterion  we  consider.  Such 
quantizer  detectors  were  presented  in  [30,3]J  under  the 
assumption  that  the  noise  distribution  is  known  and  in 
[32]  under  the  assumption  that  the  noise  distribution 
is  partially  known.  The  following  remarks  are  in  order. 

1)  With  the  proper  choice  of  detector  will 

limit  the  effect  of  large  observations  to  the  desired 
degree. 

2)  The  detector  is  insensitive  to  the  specific  shape 
of  the  distribution  since  it  does  not  discriminate 
between  the  values  of  observations  in  the  same  interval. 

3)  If  there  is  a family  of  distributions  such  that 
[F(aj)-F(aj  ] under  Hq  is  the  same  for  all  members 
of  this  family  for  j=l,...,  m,  then  the  distribution 

of  Tj^  will  be  the  same  for  any  member  of  this  family  under 
Hq  and  the  test  will  be  a constant  false  alarm 
rate  test  (CFAR)  over  this  family. 


From  the  above  remarks  and  since  we  are  concerned  with 
the  threshold  stability,  it  seems  that  the  most  appro- 
priate way  to  choose  {a^}  is  in  some  way  related  to  the 
quantiles  of  the  distribution.  This  simplifies  the 
implementation  of  the  receiver  in  the  case  of  detection 


97 


in  sea  clutter  since  the  designer  usually  has  at  least 
some  idea  about  the  quantiles  of  the  distribution.  In 
all  detectors  involved  in  this  study,  the  ^j's  were 
choosen  such  that  [F (a ^ ) -F (a i for  all  j=l,...,  m. 
To  choose  the  Ij's  three  different  methods  were  used 

a)  The  fixed  step  quantizer 

In  this  method,  1^  was  chosen  as  ^j~j* 

assume  that  n^  observations  fall  between  and 

a.,  the  test  statistic  can  be  rewritten  as 
J 


m 


T.,  = ^ y j n . 

N N -•  D 


(V-20) 


To  find  the  distribution  of  T„  we  first  find  the 

N 

distribution  of  the  n ^ ' s which  is  given  by 


Prl^i-Ni,  1=1,2,...,  . . .N_!  (V-21) 


m 


and 


r,  r.  N 

J-  Pm 

Pr[ni=Ni,  i=l,...,  m|Hj^]=N!  ) • • • 

1 * m ■ 


(V-22) 


where  P^  is  the  probability  that  one  observation  falls  in 
the  interval  j under  Using  equation  (V-20)  - (V-22)  we 


can  get  the  distribution  of  T^^, 


b)  Coherent-Optimal  quantizer; 

In  [32]  it  was  shown  that  if  the  signal  is  a dc  signal 
in  additive  noise  then  a locally  most  powerful  quantizer 
detector  of  the  above  form  can  be  constructed  using 
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1 . 
D 


^ — !/•••  / in 


(V-23) 


Although  this  is  not  the  case  here,  this  quantizer  was 
included  in  the  study  for  the  sake  of  completeness.  If 
we  assume  the  observations  are  distributed  under  Hq  as 
a lognormal  with  a=6db  then  the  values  of  the  a ^ ' s and 


1 . ' 
3 

s would 

be  given 

from 

table  (V-1) 

below 

for 

m=10 

j 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

a . 

D 

.41 

.556 

.693 

.841 

1.0 

1.10 

1.44 

1.8 

2.43 

00 

1 . 

■) 

-.613  - 

-.111 

0 

.059 

.088 

.107 

.123 

. 125 

.119 

.104 

Table 

(V-1) 

a . ' s 
3 

and 

1 . ' s for  10 
3 

level 

coherent 

optimal  quantizer  under  lognormal  distribution 


Note  that  all  the  Ij's  for  j>5  are  in  the  neighborhood 
of  .11,  and  that  the  values  of  the  Ij's  for  j^8  are 
decreasing.  This  second  remark  is  a setback  in  this 
method  since  it  implies  that  as  the  signal  increases 
may  decrease,  which  would  be  expected  to  lead  to  bad 
performance. 

c)  Semi  optimal-coherent  quantizer: 

To  correct  for  the  descending  part  in  the  l^'s,  1^ 
and  Ij^Q  were  set  equal  to  0.126  and  0.127  respectively. 
Otherwise,  this  quantizer  is  the  same  as  the  previous  one. 
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ji 

Numerical  Results 

To  evaluate  the  actual  performance  of  the  above  three 
detectors,  simulation  experiments  were  conducted  using 
lognormal  and  contaminated  normal  distributions  as  in 
Section  (V-1) . The  results  for  the  optimal  coherent 
quantizer  detector  were  as  bad  as  expected  especially 
for  low  values  of  P^,,  so  the  results  were  not  reported  here. 

Figures  (32-34)  show  the  probability  of  detection 
versus  the  signal  to  noise  ratio  (S/N)  at  different 
probabilities  of  false  alarms  P , for  the  fixed  step 

r 

quantizer  detector.  Comparing  figure  (33)  with  (34) 

we  can  see  that  the  signal  to  noise  ratio  required  to 

achieve  a certain  P,  at  the  same  P„  when  the  noise,  is 

d F 

lognormal  is  very  close  to  that  required  when  the  noise 
is  contaminated  normal. 

Figure  (35)  shows  P^  versus  (S/N)  for  the  semioptimal 
coherent  quantizer  detector.  Comparing  figure  (33)  with 
(35)  we  notice  that  the  fixed  step  quantizer  detector 
is  superior  to  the  semioptimal  coherent  one  especially 
at  higher  values  of  P^. 

Summary 

In  this  section,  three  quantizer  detectors  were 


4 


investigated.  It  was  shown  that  the  simplest  quantizer  — 
which  is  the  fixed  step  quantizer  — is  the  best  one. 

To  implement  this  quantizer  detector  we  need  only  to 


S/N  (dB) 


Probability  of  detection  vs  S/N  (dB)  for  log-normal 
distribution  (a=6dB)  and  equal  spaced  quantizer 
detector,  N=30 

Figure  32 


S/N  (dB) 

Probability  of  detection  vs  S/N  (dB)  for  log-normal 
distribution  (a=6dB)  and  equal  spaced  quantizer 
detector;  N=10 

Figure  33 


1 

,1 

;1 


S/N  (dB) 


Probability  of  detection  vs  S/N  (dB)  for  log-normal 
distribution  (o=6dB)  and  semi-optional  coherent  quantizer 

detector;  N=10 

Figure  35 
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know  some  the  distribution's  quantiles.  The  number  of 
quantiles  assumed  to  be  known  in  the  experimental  study 
here  was  9.  Since  this  number  is  large,  it  is  our 
opinion  that  some  study  should  be  devoted  to  the  estima- 
tion of  these  quantiles.  A first  step  in  this  direction 
was  taken  by  the  author  in  [33]  for  the  case  of  detection 
coherent  signals  in  symmetric  noise.  The  fixed  step 
quantizer  detector  will  be  compared  to  other  detectors 
later  in  this  chapter. 
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F(u,)  = i (V-28) 

1 n 

= n f(u^)  (V-29) 

and  n is  the  number  of  samples  from  which  the  maximum 
or  minimum  sample  is  taken.  In  [35]  Gumbel  also  gave  a 
technique  to  estimate  the  parameters  Uj^,u^,aj^  and  . 

Using  the  above  results  Milstein  et  al  [36  ] introduced 
the  "Extreme  Value  Theory  Detector"  (EVT) . The  main  idea 
in  this  detection  scheme  is  that  if  the  signal  is  fairly 
large, then  the  confusion  results  as  to  whether  a certain 
observation  is  related  to  the  lower  tail  of  the  distribu- 
tion of  signal  plus  noise  or  to  the  upper  tail  of  the 
noise  only  distribution.  To  solve  this  problem  they 
utilized  an  approximation  of  the  likelihood  ratio  depending 
on  the  distribution  approximation  given  in  equations 
(V-24)  and  (V-25) . So,  the  test  for  a single  observation 
will  be  of  the  form 


exp  [-a^{x-u„)]  .] 

A(x)  = f 3 tt  < 

exp  [-a^(x-u^)] 


(V-30) 


where  y is  the  threshold.  Now,  if  there  are  N observa- 
tions we  repeat  the  test  for  each  one  of  them  and  count 
how  many  observations  out  of  these  will  exceed  the 
threshold.  If  the  probability  of  false  alarm  for  a single 

observation  is  P_,  and  m observations  exceeded  the 
FL 

threshold,  then  the  total  P_  can  be  calculated  using  the 

r 


binomial  distribution.  This  detector  can  be  implemented 
very  simply  by  setting  two  thresholds.  The  first  is  y , 
which  would  be  exceeded  under  by  any  single  observation 
with  probability  The  second  is  K,  which  is  equal  to 

the  total  number  of  observations  exceeding  y , to  get 
the  total  probability  of  false  alarm  P . 

r 

Notice  that  this  detector  is  the  same  as  the  binary 
integrator  detector  presented  in  [29]  for  detection 
of  signals  in  log  normal  clutter.  It  was  also  presented 
in  [37] for  detection  of  signals  in  Rayleigh  noise  channels. 
In  general,  it  is  another  implementation  of  the  rank 
detector. 

A problem  which  might  appear  with  such  detectors  is 
in  the  choice  of  the  second  threshold  K.  Test  results 
show  that  the  receiver  performance  depends  on  the  choice 
of  K and  there  is  a different  optimal  choice  for  each 
distribution. 

Numerical  Results 

We  conducted  simulation  experiments  for  the  above 
detector  using  lognormal  and  contaminated  normal  distribu- 
tions as  done  before  in  the  previous  two  sections.  First 
we  considered  y to  be  constant  (it  was  changed  in  a 

very  narrow  region  to  guarantee  P_)  and  changed  the 

r 

second  threshold  K to  get  the  required  value  of  P . The 

r 

results  of  this  test  are  shown  in  figures  (36)  to  (38). 


S/N  (dB) 

Probability  of  detection  vs  S/N  (dB)  for  log-normal 
distribution  (o=6dB)  and  E.V.T.  detector,  N=10 


Figure  37 
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It  should  be  noticed  that  the  required  (S/N)  for  lognormal 

is  slightly  higher  than  that  for  contaminated  normal  to 

achieve  the  same  P , , especially  at  low  P . The  test  was 

repeated  for  lognormal  using  K=7  and  changing  y when 

N=10.  Figure  (39)  shows  the  results  of  this  test,  which 

show  improvement  in  the  receiver  performance,  especially 
-2  -3  -4 

at  P_  =10  ,10  and  10  where  K was  taken  in  the 

F 

previous  test  to  be  equal  to  4,5,6  respectively. 


Summary 

In  this  section,  the  EVT  detector  was  described  and 
it  was  shown  that  it  leads  to  the  same  detection  scheme 
as  the  binary  integration.  It  was  shown  also  that  the 
test  performance  is  sensitive  to  the  second  threshold 
K which  cannot  be  chosen  in  an (optimal  way  without 
knowing  the  exact  distribution  under  consideration. 


(V-5)  The  Trimmed-Mean  detecror ; : 

In  this  section,  we  give  a brief  description  of  the 
trimmed-mean  detectors  which  will  be  utilized  in  the  next 
section  for  comparison.  Assume  there  are  N observations, 
then  the  designer  chooses  two  numbers  n^^  and  n2  such  that 
N^n2^n^^0,  and  the  test  statistics  will  be 


T 


N 


n. 


2 

I 

i=n. 


z . 
1 


where  z,  is  the  ith  smallest  observation,  i.e.  z,<z_<. 

1 1—2— 


. <z 
— n 


2 


6 


8 


S/N  (dB) 

Probability  of  detection  vs  S/N  (dB)  for  log-normal 
distribution  and  E.V.T.  detector  (m=7) ; N=10 


Figure  39 
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One  difficulty  with  this  test  is  that  it  requires  sorting 

the  data.  The  second  difficulty  is  in  setting  the 

threshold  since  it  depends  on  the  distribution  and  there 

is  no  mathematical  expression  to  be  used.  Notice  also 
N 

that  if  nj^=n2=[^]  this  turns  out  to  be  the  median  detector 
and  if  nj^=l  and  n2=N  it  is  the  mean  detector.  On  the 
other  hand,  it  was  shown  in  [27]  that  this  detector 
performs  better  than  most  of  the  known  detectors  for 
both  lognormal  and  contaminated  normal  distributions. 

More  details  about  the  trimmed-mean  detector  can  be 
found  in  [38 ] . 

Figure  (40)  shows  the  performance  of  the  trimmed 
mean  detector  against  lognormal  [a=6db)  when  N=10,  nj^=2 
and  n2=5.  The  performance  of  this  detector  against 
contaminated  normal  is  the  same  as  the  M-detector  on 
envelope,  so  it  will  not  be  repeated  here.  It  was  found 
also  that  the  M-detector  is  better  for  both  distributions 
then  the  trimmed  mean  n^=4,  . 

(V-6)  Detector  Comparison 
a)  Performance 

Comparing  the  performance  results  of  the  chosen  detectors, 
it  should  be  noted  that  for  the  lognormal  distribution  at 
N=10  the  performance  of  the  M-detector  on  envelope,  the 
trimmed-mean  detector  (n^-2,  n2=5)  and  the  fixed  step 


S/N  (dB) 


Probability  of  detection  vs  S/N  (dB)  for  log-normal 
distribution  (a=6  dB)  and  Trimmed  mean 
detector  (nj^=2  ,n,-5)  ;N=10 

Figure  40 
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quantizer  detector  are  all  very  comparable.  In  fact 

the  graphs  for  against  (S/N)  overlap  in  most  cases, 

and  all  of  them  are  within  less  then  Idb.  Also,  the 

performance  of  the  M-detector  gets  bette  than  the 

others  as  P_  decreases  and  P,  increases.  The  E.V.T. 

F d 

detector  always  has  poorer  performance  than  the  others, 
even  for  K=7  which  was  shown  to  be  optimal  for  lognormal 
when  N=10  [29].  Figure  (41)  shows  an  example  of  these 
comparisons  at  P =10  It  shows  also  the  Chernoff 

r 

lower  bound  as  an  approximation  for  the  performance 

of  the  optimal  detector.  Both  the  M-detector  and  the 

fixed  quantizer  detectos  are  the  closest  to  this  bound, 

with  the  M-detector  better  at  higher  values  of  P^. 

At  N=30,  the  performance  of  both  the  M-detector  and 

the  fixed  step  quantizer  is  almost  the  same  with  the 

M-detector  better  at  low  P,.  The  E.V.T.  is  far  behind, 

d 

mainly  because  we  did  not  utilize  the  optimal  value  for 

K. 

For  the  contaminated  normal,  the  performance  of  the 

trimmed  mean  and  M-detector  was  the  same  as  mentioned 

before.  The  E.V.T.  was  the  worst  detectoi  in  the  group. 

The  fixed  quantizer  is  very  slightly  better  than  the 

M-detector  at  P_  = 10  but  the  M-detector  is  much 
r 

better  at  P <10 

I 


L 
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General  properties 

1)  The  performance  of  the  E.V.T.  detector  is  highly 
dependent  on  the  value  of  the  second  threshold  K and 
the  true  noise  distribution,  and  since  the  main  reason 
behind  this  study  is  the  uncertainty  about  the  noise 
distribution,  this  detector  cannot  be  recommended  for 
applications. 

2)  Although  the  trimmed-mean  detector  has  comparable 
performance  to  the  other  two  good  detectors,  implementation 
seems  impractical  since  it  requires  sorting  the  data, 

and  there  is  no  known  way  to  set  the  threshold  for  it. 

3)  The  M-detector  and  the  fixed  step  quantizer  detector 
are  the  simplest  to  implement  among  all  detectors 
studied,  and  they  also  have  the  best  perfori.  ance. 

4)  The  M-detector  has  the  advantage  over  the  fixed  step 
quantizer  detector  that  it  requires  less  information. 

It  requires  knowledge  of  only  two  quantiles,  while  the 
quantizer  detector  requires  more  quantiles. 


5)  Using  the  quantizer  detector,  P„  will  be  set  only 

r 


approximately  and  some  values  of  P„  are  unachievable 

r 


for  finite  N. 

6)  If  the  designer  is  interested  in  the  low  P^  high  P^ 
region,  the  M-detector  has  better  performance  than  the 
quantizer  detector. 


i 

J 
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It  is  our  opinion  that  the  M-detector  is  the  best 
of  all  the  detectors  involved  in  the  study,  for  the 
above  reasons.  Meanwhile,  we  believe,  as  was  mentioned 
earlier  in  this  chapter,  that  some  research  efforts 
should  be  directed  toward  finding  procedures  for  setting 
the  quantizer  parameters  that  depend  on  less  information 
than  is  used  here. 

(V-7)  Summary 

The  problem  of  detection  of  signals  with  random 
phase  in  sea  clutter-like  noise  was  treated  in  this 
chapter.  Several  detectors  were  developed  and/or 
described.  Comparing  these  detectors  as  to  performance 
and  implementation  it  was  found  that  the  M-detector 
on  envelope  is  best. 
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VI.  Performance  of  the  M-Detector  with  Dependent  Data 


In  all  the  above,  it  was  assumed  that  the  noise  samples 
were  independent.  In  this  chapter,  we  will  examine  the 
performance  of  M-detectors  in  certain  types  of  correlated  noise. 

Let  {x^}  be  a dependent  sequence  of  random  variables 
given  by 

m 

X — 5]  X . X . + w , n — 0,+l,  4-2,... 

n 1— I ^ n-i  n — — 

where  {w^}  is  a sequence  of  independent,  identically 
distributed  random  variables  with  mean  zero.  Let  the  observa- 
tions be 


y = 0 4-  X 
■^n  n 


and  let  6 satisfy 


, N 


We  want  to  show  that  6 ->^6 


a.s.  and  that  *^{6-0)  is 


asymptotically  normal. 

Our  princioal  tool  is  the  following  result  due  to 
Billingsley  [39,  P.184].  Let 


n — ^ ^ ^ 1/2,... 

n n n-1 


'^£n  " ^H^^n'^n-l'*’"'  ^n-£4-l^  ' 


v(£)  = E{|nQ-  nj^l  , £=1.2, 
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Sn  = rii+n2+-..+n„ 

where  f(-)  and  f^(')  are  measurable  functions,  is  a 

1/2 

({i-mixing  sequence  with  I <"  / and  has  zero  mean 

1/2 

and  finite  variance.  Suppose  also  that  ^[v(£)]  <"  . Then 


the  series 


= ECn?)  +2  I E(n«n^) 
U j=l  ^ J 


2 1 

converges  absolutely;  if  o >0,  then  ^ converges  in 

2 

distribution  to  N(0,o  ). 

To  apply  this  result,  we  rewrite  as 


X = y b .w  . 
n j£o  D n-3 

m 

b . = I a.  r . ^ 

3 i£i  ^ ^ 


in  iH”  1 

where  the  r^'s  are  the  roots  of  r -Aj^r  -...-A^  and  the  a^'s 
are  constants  (see  Sirvanci  and  Wolff  [40]).  Then 

“>  2 _ 


= I b a < » . 
n j=o  J ^ 


Note  that  {w^}  is  (})-mixing  since  it  is  an  independent  sequence. 


It  remains  to  show  that 


°°  2 1/2 
I .{E{l(x  )-l{x  ))  } / < oo 

£=1  ^ 


J. 

= y b .w  . 

j = 0 => 


(VI. 1) 


Lemma  (VI . 1)  ; 


If  I l(x^)-l(x2)  I £ R |xj^-X2l  for  all  Xj^,X2,  and  for 
0<R<<»,  then  (VI.  1)  holds. 


Proof : 


We  have 


1 KlXo-x^jl  - R|  I b.w..l 

D ^ 


and  thus 


;{i(x„)-i(x„,))2  < e(  b^ 


Let  a be  the  maximum  of  the  {a^}  and  r be  the  maximum  of 


the  . Then  since  | b ^ | £ mar^  , 


r , 2 ^ 2 r ,2  2j  2 2 r 

I b._<m  2.  ar  =ma  — 


21 


2 = 1 


j = l 


1-r 


and  therefore 


I {E(l(x  )-l(x  < — I = 

Z=1  ° “ /7T2  Z=1  /7^ 


Rma  ar 
w 


1-r 


1-r  (1- 


since  0<r<l. 

Note  that  if  1(*)  is  either  of  the  following,  it 
satisfies  the  requirements  of  the  above  lemma: 


a , t_>a 

l(t)  = t , -a<t<a 
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or 


l(t) 


with  Icaj  < Tr/2  . 

A 

To  prove  consistency  of  0 
of  the  above.  Define 


tan  (ca)  , t^a 
tan(ct)  , -a£t£a 
-tan(ca)  , t£-a 

we  need  a slight  generalization 


n = l(x  -0)  - E[l(x  -0)] 
n n ‘1 


Then 

2 

v(J!-)  = E{  InQ-no^l  ^ 

= E{  1 1(Xq-0)  -1(Xqj^-0)  1 ^}-{E(1(Xq-0)-1(Xqj^-0)  ) } 

< E{1(Xq-0)-1(Xq^-0)  }^1R^E(Xq-Xqj^)^  . 

and 

“ 1/2 
I [v(£)l  ^ <“  . 

1=1 

Now  from  the  above  we  have 

N 2 

_L_  y [l(x.-0)-E[l(x.-0)  ] 1 N(0,Oj^) 

/N  j=l  ^ 


00 

= E(no)  + 2 E(nonj)  • 


where 


2 

'^l 
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That  is,  for  any  e>0 


1 N 

i l(Xj-(0Q+e))  ^ EEKx^-Op+e))]  > 0 

1 N 

i I l(x.-(6Q-e))  E[l(x^-(0Q-e))]  < 0 


where  0^  satisfies 


E[l(Xi-0o)]  = 0 . 


Thus  for  any  there  is  a Nq  such  that  for  N > Nq 


% - e i *0  * = 


and  therefore  0 ->■  0^  a.s. 

Since  /N  0 _<  y if  and  only  if  — ^ 1 (x  .-y/v^)£0  , 

/n  j=l  ^ 

asymptotic  normality  follows  from  the  proof  of  Thm.  2.1.  in 
[39]  if  we  take  = l(x^-  y//N) . 

To  evaluate  the  finite  sample  size  performance  of  the 
M-detector  on  components  and  on  envelope  in  the  presence  of 
correlated  noise,  a set  of  simulations  were  conducted  using 
correlated  data.  For  the  M-detector  on  components  a correlated 
Rayleigh,  log-normal,  and  contaminated  normal  noise  were 
utilized.  To  generate  the  correlated  Rayleigh  (z)  we  utilized 
the  following  formula 
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X 


i 


0.1 


^i-1  ^ 


u . 
1 


y.  = 0. 1 y . , + V. 
■'i  -'i-l  1 

, 2 ^ 2,1/2 
2 . = (x . + y . ) 

1 1 1 


where  {u^}  and  independent  samples  from  a zero 

mean  normal  distribution  with  variance  0.99.  This  value  was 
chosen  to  make  f(x)  and  f (y)  members  of  the  class  of  p-point 
distributions  containing  the  standard  normal  distribution. 

This  will  simplify  the  comparison  with  the  independent  case 
treated  in  Chapter  2.  The  correlated  log-normal  was  generated 
as 


2^  = exp(y^) 

Yi  = 0-1  Yi_i  + 

where  {u^}  is  as  before.  The  contaminated  normal  was  generated 
as  in  the  Rayleigh  case  except  that  {u^}  and  were 

independent  scimples  from  a contaminated  normal  distribution. 

Figures  42-44  show  performance  of  the  M-detector  on 
components  (same  set  of  detectors  as  in  Section  II. 5)  for 
log-normal,  contaminated  normal,  and  Rayleigh,  respectively. 

Also  shown  is  the  performance  of  the  conventional  square-law 
detector.  The  M-detector  is  still  superior  to  the  square-law 
detector  for  both  log-normal  and  contaminated  normal  cases, 
and  it  is  less  than  IdB  worse  than  the  square- law  detector  for 
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the  Rayleigh  distribution.  Comparing  these  graphs  with 
Figs.  2-4,  it  can  be  seen  that  there  is  a loss  of  power  that 
achieves  a maximum  of  about  IdB  at  P„  = 10  ^ when  tl:e  noise 
is  Rayleigh.  The  minimum  power  loss  occurs  in  the  log-normal 
case. 

The  graphs  do  not  tell  us  how  much  power  is  lost  if  the 
threshold  is  set  assuming  correlated  noise  when  in  fact  the 
noise  is  independent.  This  power  loss  was  calculated  and 
found  to  be  less  than  1.5dB,  the  maximum  loss  occuring  in 
the  Rayleigh  case. 

To  test  the  performance  of  the  M-detector  on  envelope 
for  correlated  noise  the  log-normal  and  contaminated  normal 
cases  were  considered.  Figs.  45-46  show  the  probability  of 
detection  vs.  signal-to-noise  ratio  (S/N)  for  these  cases 
at  N=10.  Also  shown  for  comparison  are  the  corresponding 
results  for  the  independent  case.  The  loss  of  power  is  always 
less  than  IdB  in  both  cases,  increasing  with  decreasing  P„. 

.C 

If  the  threshold  was  set  to  give  P„  = 10  ^ assuming  correlated 

r 

noise,  then  the  actual  P would  be  10  ' if  in  fact  the  noise 

r 

was  uncorrelated.  This  is  equivalent  to  a IdB  power  loss. 

For  larger  values  of  P^,  the  loss  is  smaller. 

Finally,  one  may  do  better  by  choosing  an  !(•)  other 
than  the  tangent  function.  However,  an  intelligent  choice 
would  require  knowledge  of  the  correlation,  which  is  usually 
unknown  a priori. 
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Figure  45 
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on  envelope  and  correlated  contaminated  normal 
noise  ( — ) and  independent  contaminated  normal  (- 
N=10. 
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VI.  Summary  and  Conclusion 

In  this  paper  the  problem  of  robust  detection  of 
non-coherent  signals  in  noise  is  solved  under  the 
assumption  that  the  noise  distribution  is  unknown,  but 
a member  of  a known  class  of  distribution  functions. 

This  problem  was  divided  into  three  main  categories: 

a)  Partially-coherent  signals 

b)  Unknown  frequency  and  phase 

c)  Random  phase 

- In  the  partially-coherent  case,  it  was  assumed  that 

the  unknown  phase  is  constant  over  each  observation 
period  but  changes  randomly  from  one  observation  period 
to  another.  We  have  introduced  two  general  methods  for 
the  design  of  detectors  which  are  asymptotically  most 
robust,  in  the  sense  of  guaranteeing  a non-trivial  lower 
bound  on  the  receiver  performance  over  the  class  of 
distributions  under  consideration.  It  is  shown  that 
the  design  of  any  of  these  detectors  requires  the 
existence  of  a density  function  of  minimum  Fisher 
information  number  in  the  class  of  distributions.  Two 
detectors  were  designed  for  the  special  case  in  which 
the  distribution  of  the  noise  inphase  and  quadrature 
phase  components  are  members  of  the  class  of  p-point 
I distributions.  Simulation  results  for  finite  sample 

sizes  were  also  given  for  different  distributions.  It 
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was  found  that  the  M-detector  performs  bettor  than  both 
A-detector  and  the  square  law  detector.  The  above 
results  were  extended  also  to  the  second  case  of  unknown 
frequency  and  phase. 

In  the  third  case,  it  was  assumed  that  the  signal 
phase  changes  randomly  from  one  observation  to  another. 

Three  detectors  were  presented,  the  M-detector  on 
envelope,  the  quantizer  detector,  and  the  extreme  value 
theory  (E.V.T.)  detector.  It  was  found  that  the  M-detector 
is  the  best  among  the  above  three  and  the  trimmed  mean 
detector,  judging  from  performance  and  ease  of  implementation 
in  the  presence  of  limited  information  about  the  noise 
distribution.  All  of  the  above  detectors  assume 
knowledge  of  some  quantiles.  The  quantizer  detector  is 
the  second  best. 

In  general,  we  recommend  the  M-detector  on  component 
or  on  envelope  when  the  observation  distribution  is  not 
completely  known,  because  of  their  threshold  stability 
property  and  because  they  protect  the  decision  rule 
against  bad  observations. 
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